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Abstract. We develop the approach to the problem of integrable discretization based on the 
notion of r-matrix hierarchies. One of its basic features is the coincidence of Lax matrices of 
discretized systems with the Lax matrices of the underlying continuous time systems. A common 
feature of the discretizations obtained in this approach is non-locality. We demonstrate how to 
overcome this drawback. Namely, we introduce the notion of localizing changes of variables and 
construct such changes of variables for a large number of examples, including the Toda and the 
relativistic Toda lattices, the Volterra lattice and its integrable perturbation, the second flows 
of the Toda and of the Volterra hierarchies, the modified Volterra lattice, the Belov-Chaltikian 
lattice, the Bogoyavlensky lattices, the Bruschi-Ragnisco lattice. We also introduce a novel class 
of constrained lattice KP systems, discretize all of them, and find the corresponding localizing 
change of variables. Pulling back the differential equations of motion under the localizing changes 
of variables, we find also (sometimes novel) integrable one-parameter deformations of integrable 
lattice systems. Poisson properties of the localizing changes of variables are also studied: they 
produce interesting one-parameter deformations of the known Poisson algebras. 
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1 Introduction 

This paper deals with some aspects of the following general problem: how to discretize one 
or several of independent variables in a given integrable system, maintaining the integrability 
property? We call this the problem of integrable discretization. 

To assure the coincidence of the qualitative properties of the discretized models with that 
of the continuous ones becomes one of the central ideas of the modern numerical analysis, 
which therefore comes to a close interplay with different aspects of the theory of dynamical 
systems. One of the most advanced examples of this approach is the symplectic integration, 
on which recently the first monograph appeared ||SSC| . 



The problem of integrable discretization constitutes another aspect of this general line 
of thinking. It arose in the course of development of the theory of solitons. This theory, 
born exactly 30 years ago, has grown so tremendously, that it is difficult to keep an overview 
of the whole variety of different notions and contexts of integrability, not to say about the 
concrete results. (Recall that recently there appeared a thick book under the title "What is 
integrability?" ||Zakh|| ). 

Correspondingly, various approaches to the problem of integrable discretization are cur- 
rently available. They began to be discussed sporadically in the soliton literature starting 
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from the mid-70s. Following ones should be mentioned: 



2. 



4. 



An approach based on the representation of an integrable system as a compatibility 
condition of two auxiliary linear problems. A natural proposition is to discretize one 
or the both of them ||AL1|] . This, however, can be made in a great variety of ways, 
cf., for instance, different spatial discretizations of the nonlinear Schrodinger equation 
and of the sine-Gordon equation, found in [|AL1|| and in |]IK| . One attempt of fixing 



discretization crystallized with the development of the Hamiltonian approach. Namely, 
Faddeev and Takhtajan, based on the experience of the Leningrad soliton school, for- 
mulated in |[l'"r2|| the following rule for a transition from models with one continuous 
space variable to lattice models: the r-matrix should be preserved, the linear Poisson 
bracket being replaced by the quadratic one. See |p:'"l'2[ | for a collection of examples 
showing productivity of this approach. 

One of the most intriguing and universal approaches is the Hirota's one 0], based on 
the notion of the r-function and on a bilinear representation of integrable systems. 
It seems to be able to produce discrete versions of the majority of soliton equations, 
but still remains somewhat mysterious, and the mechanism behind it is yet to be 
fully understood. One successful way to do it was proposed in IPJM]] , where also a 



large number of integrable discretizations was derived. Among the most interesting 
products of this approach is the so called Hirota-Miwa equation [[Hl|| , \fA^ , which is 
sometimes claimed to contain "everything" , i.e. the majority if not all soliton equations 



(continuous and discrete) are particular or limiting cases of this single equation, cf. 



mm 



A fruitful method is based on the "direct linearization" 
[[NPCQII , [[NQ| . Its basic idea is to derive integrable nonlinear differential equations 



which are satisfied by the solutions of certain linear integral equations. A large variety 
of continuous and discrete soliton equations has been obtained on this way. 



Approach based on the variational principle (discrete Lagrangian equations 
with matrix factorizations 



^ combined 

|MV|] , |PLT2|| . Historically, it was the work of Veselov 



and Moser that consolidated the more or less isolated results to a separate branch of 
the theory of integrable systems. 

Considering stationary and restricted flows of soliton hierarchies, closely related to the 
"nonlinearization" of spectral problems, often leads to interesting discrete equations 



6. Differential equations describing various geometric problem (surfaces of the constant 
mean curvature, motion of the curve in the space, etc.) turn out to be integrable |pym|| . 
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B2|| . Correspondingly, a discretization of geometric notions naturally leads to 



discrete integrable equations | |BPl| j, [ |BP2| ] , [ pS| | , |^ 



7. There exist integrable discretizations which belong to the most beautiful examples, but 
were derived by guess, without any systematical approach [ pT| , 



8. Last but not least we menton an approach to the temporal discretization in which the 
auxiliary spectral problem is not discretized at all. In other words, the basic feature 
of this approach is maintaining the Lax matrix of the continuous time system. The 



first example of this approach is the work by Ablowitz-Ladik ||AL2|| , futher developed 
towards the practical algorithms in [[TA|| . This feature was also put in the basis of 
the work by Gibbons and Kupershmidt ||GK]| , [|K2|| . The discretizations found in all 
these papers were somewhat unsatisfactory from the esthetical point of view, namely 
they suffered from being nonlocal, as opposed to the underlying continuous time sys- 
tems. Morever, these authors did not recognize the connection with the factorization 
problem, which did not allow them to identify their discretizations as certain members 
of the corresponding hierarchies and to establish the Poisson properties of these dis- 
cretizations. This led Gibbons and Kupershmidt to call this method "the method of 
the bizarre ansatz". 

Recently the author pushed forward the last mentioned approach to the problem of 
integrable discretization, putting it in a connection with the r-matrix theory of integrable 
hierarchies (see ||RSTS|| for a review of this theory). In this context the method could be 
understood properly, and became rather natural and simple. It was applied to a number of 
integrable lattice systems ||56|- ||S11[] , ||S13|| , ||S14|| . Its clear advantage is universality. The 
method is in principle applicable to any system admitting an r-matrix interpretation, which 
is the common feature of the great majority of the known integrable systems. As for the 
drawback of nonlocality, there exist several ways to repair it. The first one, connected with 
the notion of discrete time Newtonian equations of motion, was followed in KB], P^, ||S10 



| S11|| . A splitting of complicated flows into superpositions of more simple ones was used 



m 



jl3|| , | [S14|| . The present work is devoted to another way connected with the so called 
localizing changes of variables. 

The paper has the following structure. In Sect. |^ we give an accurate formulation of 
the problem of integrable discretization. Sect. |^,^ are devoted to a general framework 
of integrable r-matrix hierarchies of Lax equations on associative algebras. In Sect. ^ we 
formulate a general recipe of integrable discretization, and in Sect, ^we introduce the notion 
of localizing changes of variables and discuss their general properties. Sect. |^ contains the 
description of algebras used for the analysis of all integrable lattice systems in this paper. 
The rest of the paper is devoted to a detailed elaboration of a number of examples, including 
the most prominent ones, such as the Toda lattice and the Volterra lattice, and less well 
known ones, such as the Belov-Chaltikian lattice. 
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2 The problem of integrable discretization 

Let us formulate the problem of integrable discretization more precisely. Let X he a Poisson 
manifold with a Poisson bracket {■,■}. Let if be a completely integrable Hamilton function 
on X, i.e. the system 

x = {H,x} (2.1) 

possesses many enough functionally independent integrals Ik{x) in involution. 
The problem consists in finding a map X ^ X described by a formula 

x = ^{x;h) (2.2) 

depending on a small parameter /i > 0, and satisfying the following requirements: 

1. The map ( |2.2| ) is a discrete time approximation for the flow ( p^.lD in the following sense: 

^{x;K)=x^h{R,x}^0{h^) (2.3) 

(of course, one might require also a higher order of approximation). In all our consid- 
erations and formulas we pay a special attention to a simple and transparent control 
of the continuous limit — > 0. 



2. The map ( |2.2D is Poisson with respect to the bracket {■,■} on X or with respect to 
some its deformation {■, -jh such that {■, - jh = {■,■} + 0(/i). 



3. The map ( p.2|) is integrable, i.e. possesses the necessary number of independent 
integrals in involution Ik{x]h) approximating the integrals of the original system: 
Ik{x;h) = lk{x) + 0{h). 

3 Lax representations 

Our approach to the problem of integrable discretization is applicable to any system allowing 
an r-matrix interpretation, but we formulate the basic recipe in a simplified form, applicable 
to systems with a Lax representation of one of the following types: 

L= [L,7r+(/(L))] =-[L,7r_(/(L))] (3.1) 

or 

4- = L,.-7r+(/(T,_i))-vr+(/(T,))-L, = - L, ■ vr. (/(T,_i)) + vr_ (/(T,)) ■ (3.2) 
Let us discuss the notations. 
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Let g be an associative algebra, supplied with a nondegenerate scalar product, which 
allows to identify g* with g. One can introduce in g the structure of Lie algebra in a 
standard way. Let g^^ g_ be two subalgebras such that as a vector space gr is a direct sum 
9 = 9+® 9 ■ Denote hy ti± : g ^ g^ the corresponding projections. Finally, let f : g ^ g 



be an Ad*-covariant function on gi, and let L stand for a generic element of g. Then (|0 
is a certain differential equation on g. 

Further, let g = ®Y=i fl* be a direct product of m copies of the algebra g. A generic 
element of g is denoted hy L = (Li, . . . , L^)- We use also the notation 

Tj = Tj{L) = Lj ■ . . . ■ Li ■ Lm • . . . • Lj^i (3.3) 

Then (|3.2| ) is a certain differential equation on g. Such equations are sometimes called Lax 
triads. 

One says that (|3.1|) , resp. (|3.2|), is a Lax representation of the flow if there exists a 

map L : X g (resp. L : X g) such that the former equations of motion are equivalent 
to the latter ones. Let us stress that when considering equations ( |3.1| ), resp. ( p.2| ) in the 
role of Lax representation, the letter L (resp. L) does not stand for a generic element of 
the corresponding algebra any more; rather, it represents the elements of the images of the 
maps L : X g and L : X g, correspondingly. The elements L{x), resp. L{x) (and the 
map L, resp. L, itself) are called Lax matrices. 

Equations ( |3.1| ) and (|3.2|) have several remarkable features. First of all, they are Hamil- 
tonian with respect to a certain linear r-matrix Poisson bracket on g (resp. on g) 



Moreover, under some natural conditions (for example, if gj^ and g_ serve as orthogonal 
complements to each other) one can define quadratic and cubic r-matrix brackets on gi, com- 
patible with the linear one, such that the equations above are Hamiltonian with respect to 
all three brackets |STS|, |Lg, |0l], [|D, W^ . 



The Lax representations are especially useful, if the maps L, resp. -L, corresponding to 
the Lax matrices, are Poisson with respect to one of the above mentioned r-matrix brackets; 
then the manifold consisting of the Lax matrices is a Poisson submanifold. In such cases one 
says that the Lax representation admits an r-matrix interpretation. 



4 Factorization theorems 



As a further remarkable feature of the equations and ( p.2|) we consider the possibility 
to solve them explicitly in terms of a certain factorization problem in the Lie group G 
corresponding to g [ ^'LS|| , [[RS'LS|| . (Actually, this can be done even in a more general 
situation of hierarchies governed by /^-operators satisfying the so-called modified Yang- 
Baxter equation, see |[RSTS|| ). The factorization problem is described by the equation 

n+(f/)n_(f/), f/GG, n±(t/)GG± (4.1) 



U 
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where G± are two subgroups of G with the Lie algebras gi^, respectively. This problem 
has a unique solution in a certain neighbourhood of the group unit. In what follows we 
suppose that G is a matrix group, and write the coadjoint action of the group elements on 
g* ~ gf as a conjugation by the corresponding matrices. [In this context we write Tl^{U) 

for {T{±{U)^ ]. Correspondingly, we call Ad*-covariant functions g ^ g also "conjugation 
covariant" . This notation has an additional advantage of being applicable also to functions 
g^G. 

For the history of the following fundamental theorem and its different proofs the reader 
is referred to [[RSTSf . 

Theorem 4.1 Let f ■ g ^ g be a conjugation covariant function. Then the solution of the 
differential equation with the initial condition L(0) = Lq is given, at least for t small 
enough, by 

Lit) = n;i(^e^/(^o)^ Lon+(^e^/(^o)^ = n_(^e^/(^o)^ Lon:i(^e^/(^o)j (4.2) 
Proof. We give a proof based on direct and simple verification. Denote 

so that 

gt/(Lo) = c{t)n{t), C{t) e G+, n{t) e G_ (4.3) 

Now we set 

L{t) = C-\t) Lq C{t) = TZ{t) Lo TZ'\t) (4.4) 

(these two expressions for L(t) are equal due to Ad*-covariance of f{L)), and check by direct 
calculation that this L{t) satisfies the differential equation ( p.l|) . The theorem will follow by 
the uniqueness of solution. We see immediately that L(t) satisfies the following Lax type 
equation: 

L = [L,c~^c] = -[L,nn-^] 

and it remains to show that 

C-'C = 7r+{f{L)), 7^7^-l = 7^_(/(L)) 

Since, obviously, C^^C G g^, IZIZ^^ G g_, we need to demonstrate only that 

C-^C + nn-^ = f{L) (4.5) 

To do this, we differentiate ( [4.3|) and derive, using Ad*-covariance of / and the definition 
(13): 

Cn + Cn = e^/(^o)/(Lo) = C 7^/(Lo) =Cf{L)n 
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This is equivalent to ( ^l5|) . ■ 

For an arbitrary conjugation covariant function F : g ^ G one can define tfie map 
■ g ^ g according to the formula 

L = Bf{L) = ni^ ■ L ■ n+ (F(L)) = n_ (F(L)) ■ L ■ n:^ (F(L)) (4.6) 



("B" is for "Backlund"). Theorem [4.1| shows that the flows defined by the differential 

equations consist of maps having such a form with F[L) = . The very remarkable 

feature of such maps is their commutativity for different F's. 

Theorem 4.2 For two arbitrary conjugation covariant functions Fi,F2 : g ^ G 

B^.oB^, =B^,p, (4.7) 

and therefore the maps Bi?^, Bp^ commute. 

Proof. Denote: 

Li = Bf,{L), L2 = Bp, o B^,(L) = Bp,{Li) 
So, by definition we have: 

Li = C^^LCi = 7^lL7^^^ L2 = L2^UL2 = U^LiTI^^ (4.8) 

where the matrices Ci G G+, IZi G G_ {i = 1, 2) come from the following factorizations: 

F^{L)=C{Ri, F2(Li) =£27^2 

From (|4.8| ) we have: 

L2 = C-^LC = nLn-\ where £ = C1C2 E G+, TZ = 7^27^l G G_ (4.9) 
Now the following chain of equalities holds: 

cn = £l£27^27^l = c^F2{L^)n^ = c^F2{c^^LC^)ni = F2{l)C{r^ = F2{l)f^{l) 

In view of ( [4. 91 ) we get: 

C = n+(F2(L)Fi(L)), n = n_(F2(L)Fi(L)) 

and the theorem is proved. ■ 

Theorem [4.2| implies that the flows of two arbitrary differential equations of the form 
( ^.1| ) commute. Another important consequence of Theorem is the following discrete- 



time counterpart of Theorem f4.1|, going back to [By 



Theorem 4.3 Let F : g G he a conjugation covariant function. Then the solution of the 
difference equation 

L = U-\F{L)) ■ L ■ n+(F(L)) = n_(F(L)) ■ L ■ Uz\F{L)) (4.10) 

where L = L{n), L = L{n + 1), with the initial condition L{0) = Lq, is given by 

L{n) = n;i (F"(Lo)) ■ Lo ■ n+ (F"(Lo)) = n_ (F^iLo)) ■ L, ■ Uz' (F"(Lo)) (4.11) 

Proof. From Theorem [4.2| there foUows by induction that (B^)" = B^n. ■ 

Comparing the formulas ([4.11|) , ( [4.2|) , we see that the map ( [4.10|) is the time h shift along 
the trajectories of the flow (|3.1|) with 

f{L) = h-Hog{F{L)) 

The above results are purely kinematic, in the sense that no additional Hamiltonian 
structure is necessary neither to formulate nor to prove them. However, as mentioned above, 
the equations ( p.lj) often admit a Hamiltonian or even a multi-Hamiltonian interpretation. 
If this is the case, then we get some useful additional information. In particular, all maps 
( [4.10| ) are Poisson with respect to the invariant Poisson bracket of the hierarchy ( p.l|) , being 
shifts along the trajectories of Hamiltonian flows. Further, if the set of Lax matrices L{X) 
for the system at hand forms a Poisson submanifold for one of the r-matrix brackets on g, 
then this manifold is left invariant by the flows (|3.1|) and by the maps ( [4.10| ). The functions 
on X of the form I o L, where / are conjugation invariants of g, are integrals of motion of 
the corresponding systems and in involution with respect to {■,■}. 

We close this section by giving analogous results for the Lax equations on the direct 
products g = <S>]Li 9- 

Theorem 4.4 For a conjugation covariant function f ■ g ^ g the solution of the equation 
( p.2|) with the initial value L{0) is given, at least for t small enough, by the formula 



L,{t) = H;1 (e^/(^.(0)) j . L,(0) ■ H^ (e^/^-iW) 

= H. fe^/(^.(0))) . L,(0) ■ H:^ fe^/(^.-i(0))) (4.12) 



Theorem 4.5 For a conjugation covariant function F : g ^ G consider the following 
system of difference equations on g.- 

L, = H;1(F(T,)) ■ L, . H+(F(T,_i)) = n_(F(T,)) ■ L, ■ Uz\F{T,^,)) (4.13) 

Its solution with the initial value L{0) is given by the formula 

LM) = h;i(f"(t,(o)))-l,(o)-h+(f"(t,_i(o))) 

= H_(F"(T,(0))) ■ L,{0) ■ H:1(F"(T,_i(0))) (4.14) 
The proofs are kinematic and absolutely parallel to the case of the "small" algebras g. 
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5 Recipe for integrable discretization 



The results of the previous section inspire the following recipe for integrable discretization, 



clearly formulated for the first time in ||S6|| , B7| 



Recipe. Suppose you are looking for an integrable discretization of an integrable system 
yil) allowing a Lax representation of the form (|3.1|). Then as a solution of your task you 



may take the difference equation ( |4.1CI|) with the same Lax matrix L and some conjugation 



covariant function F : g ^ G such that 

F{L) =1 + hf{L)+0{h^) 

Analogously, if your system has a Lax representation of the form ( ^.21 ) on the algebra g, 
then you may take as its integrable discretization the difference Lax equation ([4.13|) with F 
as above. 

Of course, this prescription makes sense only if the corresponding factors n±(F(L)) [resp. 
n±(F(Tj))] admit more or less explicit expressions, allowing to write down the corresponding 
difference equations in a more or less closed form. The choice of F is a transcendent problem, 
which however turns out to be solvable for many of (hopefully, for the majority of or even 
for all) the known integrable systems. The simplest possible choice F{L) = / + hf{L) works 
perfectly well for a vast set of examples considered below. 

Let us stress the advantages of this approach to the problem of integrable discretization. 

• Although we formulate our recipe only for systems with Lax representations of the 
particular form, it is in fact more universal. Almost without changes it may be applied 
to any system whose Lax representation is governed by an /^-operator satisfying the 
modified Yang-Baxter equation. 

• The discretizations obtained in this way share the Lax matrix and therefore the inte- 
grals of motion with their underlying continuous time systems. 

• If the Lax representation [resp. ( p.2| )] allows an r-matrix interpretation, then our 
discretizations share also the invariant Poisson bracket with the underlying continuous 
time systems. In particular, if the Lax matrices L [resp. L] form a Poisson submanifold 
for some r-matrix bracket, then this submanifold is left invariant by the corresponding 
Poisson map ( [4.101 ) [resp. ( |4.13| )]. 

• The initial value problem for our discrete time equations can be solved in terms of the 
same factorization in a Lie group as the initial value problem for the continuous time 
system. 

• Interpolating Hamiltonian flows also belong to the set of granted by-products of this 
approach. 
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6 Localizing changes of variables 



Along with these advantageous properties our recipe has also an important drawback: it 
produces, as a rule, nonlocal difference equations, when applied to lattice systems with 
local interactions. Under locality we understand the following property: in some coordinates 
(xi, . . . ,Xn) on X the equations of motion (p.l|) have the form 



ik = 4>k{x) = 4>k{xk, Xk±i, Xk±s) (6.1) 

with a fixed s G N. Nonlocal difference equations produced by our scheme have the form 

Xk = Xk + h^k{x;h), ^k{x,0) = (j)k{x) (6.2) 

where $fc depends explicitly on all Xj, not only on 2s nearest neighbours of Xk- The aim of 
the present paper is to demonstrate on a large number of examples how this drawback can 
be overcome, i.e. how to bring the latter difference equations into a local form; the price we 
have to pay is that they become implicit. 

The general strategy will be to find localizing changes of variables A'(x) i-^ X{x) such 
that in the variables x the map (|6.2|) may be written as 



Xfc = Xfc + /i^fe(x,x; /i), '^k{x,x;0) = 4>k{x) (6.3) 

where depends only on the x^'s and Xj's with correct indices \j — k\ < s. Such implicit 
local equations of motion are much better suited for the purposes of numerical simulation 
and are much more satisfactory from the esthetical point of view. Moreover, in all our 
examples the functions actually depend only on x^'s with k < j < k + s and on x^-'s with 
k — s < j < k, which makes the practical implementation of the corresponding difference 
equations even more effective (if, for instance, one uses the Newton's iterative method to solve 
(O) for X, then one has to solve only linear systems whose matrices are triangular and have 



a band structure, i.e. only s nonzero diagonals). Further, it has to be remarked that, when 
considered as equations on the lattice (t, k), the equations ( |6.3| ) often allow transformations 
of independent variables (mixing t and k) bringing these equations into local explicit form 
(cf. [pNQ] ). This last remark will be the subject of a separate publication. 



It is by no means evident that such localizing changes of variables exist, but we give in 
this paper a large number of examples which, hopefully, will convince the reader that this is 
indeed the case. The probably first examples appeared in the context of the Bogoyavlensky 
lattices in ||S^| , and will be reproduced here in a clarified form for the sake of completeness. 

The localizing changes of variables turn out to have many additional remarkable proper- 
ties. They are always given by the formulas 

Xfe = Xfc + hEk{yi; h) (6.4) 
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with local functions S^. However, the inverse change of variables is always nonlocal. There- 
fore nothing guarantees a priori that the pull-back of the differential equations of motion 



( |6.1|) under the change of variables ( |6.4|) will be given by local formulas. Nevertheless, this 
turns out to be the case. This gives a way of producing (sometimes novel) one-parameter 
families of integrable local deformations of lattice systems (see ||K1|| for a general concept 
and some examples of integrable deformations). 

The system ( |6.1| ) often admits one or several invariant local Poisson brackets. Nothing 
guarantees a priori that the pull-backs of these brackets under the change of variables (|6.4| ) 
are also given by local formulas. Indeed, as a rule these pull-backs are non-local. However, in 
the multi-Hamiltonian case it often turns out that pull-backs of certain linear combinations 
of invariant Poisson brackets are local again! 

These facts still wait to be completely understood. It seems that the remarkable prop- 
erties of the localizing maps have the same nature as that of the Miura maps ("miraculous 
cancellations"). Moreover, actually our localizing changes of variables are Miura maps, and 
in this image some of them already appeared in ||K1|] . However, the observation that they 
bring integrable discretizations into a local form, seems to be completely new. The scope of 
the present work is restricted to elaborating a large number of examples of localizing changes 
of variables, along with their Poisson properties, in a hope to attract the attention of the 
soliton community to these fascinating and beautiful objects. 



7 Basic algebras and operators 

Two concrete algebras play the basic role in our presentation. They are well suited to de- 
scribe various lattice systems with the so called open-end and periodic boundary conditions, 
respectively. Here are the relevant definitions. 

For the open-end case we always set g = gl{N), the algebra of N x N matrices with the 
usual matrix product, the Lie bracket [u,v] = uv — vu, and the nondegenerate bi-invariant 
scalar product {u, v) = tT{u ■ v). As a linear space, g may be presented as a direct sum 

g = 9+®g 

where the subalgebras g_^_ and g_ consist of lower triangular and of strictly upper triangular 
matrices, respectively. 

The Lie group G corresponding to the Lie algebra g is GL{N), the group of N x N 
nondegenerate matrices. The subgroups G+, G_ corresponding to the Lie algebras g^, g_ 
consist of nondegenerate lower triangular matrices and of upper triangular matrices with unit 
diagonal, respectively. The H_|_H_ factorization is well known in the linear algebra under the 
name of the LU factorization. 

In the periodic case we always choose as gr a certain twisted loop algebra over gl{N). A 
loop algebra over gl{N) is an algebra of Laurent polynomials with coefficients from gl{N) 



13 



and a natural commutator [u\^,v\''] = [u,v]X^~^''. Our twisted algebra fli is a subalgebra 
singled out by the additional condition 

g = {u{X) e ^/(A^)[A, A"^] : flu{X)n-^ = u{uX)} , 
where fl = diag(l,a;, . . . ,u!^~^), cu = exp{27Ti/N). In other words, elements of g satisfy 

^(A)=E E ^M?^.^ (7-1) 

p j — k=p 
(modJV) 

(Here and below Ej^ stands for the matrix whose only nonzero entry is on the intersection 
of the jth row and the kth column and is equal to 1). The nondegenerate bi-invariant scalar 
product is chosen as 

MA),^;(A))=trKA)-^;(A))o (7.2) 

the subscript denoting the free term of the formal Laurent series. This scalar product 
allows to identify g* with g. 

As a linear space, g is again a direct sum 

with the subalgebras 

»+ = 0A'g„ 9- = ®X'g, (7.3) 

fe>0 fe<0 

The group G corresponding to the Lie algebra gr is a twisted loop group, consisting of 
G'L(A^)-valued functions U{X) of the complex parameter A, regular in CP^\{0, oo} and 
satisfying QU{X)Q^^ = U{ujX). Its subgroups and G_ corresponding to the Lie algebras 
gr_,_ and g_, are singled out by the following conditions: 

• [/(A) e G+ are regular in the neighbourhood of A = 0; 

• C/(A) e G- are regular in the neighbourhood of A = oo and U{oo) — I. 

We call the corresponding n+n_ factorization the generalized LU factorization. It is uniquely 
defined in a certain neighbourhood of the unit element of G. As opposed to the open-end 
case, finding the generalized LU factorization is a problem of the Riemann-Hilbert type 
which is solved in terms of algebraic geometry rather than in terms of linear algebra. 



14 



8 Toda lattice 



8.1 Equations of motion and tri— Hamiltonian structure 

The equations of motion of the Toda lattice (hereafter TL) read: 

Ofc = afc(6fc+i - bk = ak-ak-i, I < k < N, (8.1) 

with one of the two types of boundary conditions: open-end (oq = oa? = 0), or periodic (all 
subscripts are taken (mod N), so that = a^, Bn+i = bi). 

The phase space of the TL in the case of the periodic boundary conditions is 

r = R^''{h,ai,...,bN,aN) (8.2) 

There exist three compatible local Poisson brackets on T such that the system TL is Hamil- 
tonian with respect to each one of them ||K1|| , see also Q. We adopt once and forever 



the following conventions: the Poisson brackets will be defined by writing down all nonvan- 
ishing brackets between the coordinate functions; the indices in the corresponding formulas 
are taken (mod A^). 

The "linear" Poisson structure on T is defined by the brackets 

{bk,(ik}i = —o^k, {flfej &fe+i}i = — Ofc (8.3) 
the corresponding Hamilton function for the flow TL is given by: 

1 Af N 

H2(a,6) = -E&'. + E«^ (8-4) 

^ k=i k=i 

The "quadratic" Poisson structure has the following definition: 

{bk,0,k}2 = —(^kbk, {0'k,bk+l}2 = —o^kbk+i 



{Ofc,afc+l}2 — —dk+i^k, {bk,bk+i}2 — —dk 

The Hamilton function generating TL in this bracket is: 



1.5) 



N 

Hi(a,6) = ^6fc (8.6) 

k=l 

Finally, the "cubic" bracket on T is given by the relations 

{bk,ak}3 = -ak{bl + ttk), {ak,bk+i}3 = -ak{bl^^ + ak), 

{afe,ctfe+i}3 = —'^ttkCik+ibk+i, {bk,bk+i}3 = —o,k{bk + bk+i), (8.7) 

Wk,bk+2}3 = —O-kO-k+l, {bk,0'k+l]3 = —O-kO-k+l 
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The expression for the corresponding Hamilton function, suitable in both the periodic and 
the open-end case, is nonlocal in the coordinates (a, 6). However, in the periodic case one 
has an alternative Hamilton function: 

1 ^ 

Ho(a,&) = -ElogK) (8.8) 
^ k=i 



8.2 Lax representation 

The Lax representation of the Toda lattice |^, |Q lives in the algebra g introduced in Sect. 
0. Actually, there exist different versions of the Lax representation connected with different 
ways to represent the algebra g a.s a. direct sum of its two subalgebras [pLTl|| . We discuss 
here only the LU version which corresponds in the open-end case to the decomposition of an 
arbitrary matrix into the sum of a lower triangular and a strictly upper triangular matrices. 
The Lax matrix T : T \—>- g of the TL corresponding to the generalized LU decomposition 

is: 

N N N 

T{a, b, A) = A"^ J2 (^kEk,k+i + J2 bkEk,k + ^Yl ^k+i,k (8.9) 

fc=l k=l k=l 

We use here and below a convention according to which in the periodic case i?7v+i,Af = Ei^n, 
En,n+i = En,u in the open-end case En+i,n = En,n+i = and we may set A = L 
The equations of motion ( ^.11) are equivalent to the Lax equations 



f = [T,B+] = -[T,B^] (8.10) 



with 



N N 

5+(a,6,A) = 7r+(T(a,6,A)) = J2bkEk,k + Ek+i,k (8.11) 

k=l k=l 

N 

B_{a,b,X) = 7r_(T(a,6,A)) = \-' ^ ('kEk,k+i (8.12) 

fc=i 

where it± : g ^ g± are the projections to the subalgebras g^ defined as in Sect. 1^ (the 
generalized LU decomposition). 

Spectral invariants of the Lax matrix T{a, b, A) serve as integrals of motion of this sys- 
tem. Note that all Hamilton functions in different Hamiltonian formulations belong to these 
spectral invariants. For instance, 

H2(a,6) = ^{tTT\a,b,X))^, Hi(a,6) = (trT(a, 6, A))^ 
where the subscript " 0" is used to denote the free term of the corresponding Laurent series. 
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All spectral invariants turn out to be in involution with respect to each of the Poisson 
brackets ( ^.31 ), ( |8.5| ), ( ^.T] ). Most directly it follows from the r-matrix interpretation of the 
Lax equation (|8.10|) , which can be given for all three brackets ||AM|| , |PLT1|| , ||0R|| , [^ , 



MP 



8.3 Discretization 

In order to find an integrable time discretization for the flow TL, we apply the recipe of Sect. 
^ with F{T) = I + hT, i.e. we take as a solution of this problem the map described by the 
discrete time Lax equation 

f = B-^TB+ = B_TBZ^ with B± = Il±{I + hT) (8.13) 



Theorem 8.1 ||S6|| (see also ||GK|| ). The discrete time Lax equation ( |8.13| ) is equivalent to 
the map (a, b) ^— (a, b) described by the following equations: 

ak = ak^^, h = bk + h[^ -^^^] (8.14) 

Pk \Pk Pk-1 / 

where the functions (3k = Pkia, b) = 1 + 0{h) are uniquely defined by the recurrent relation 

I3, = l + hbk~^^-p^ (8.15) 

Pk-l 

and have the asymptotics 

(3k = l + hbk + 0{h^) (8.16) 
Remark. The matrices B± have the following expressions: 

N N 

B+{a, b, A) = n+(/ + HT) = PkEk,k + hXj2 Ek+i,k (8.17) 

k=l k=l 
N 

B.{a, 6, A) = n_ (/ + /iT) = / + h\-^ J2 -^Ek,k+i (8.18) 

k=i 

Proof. The bi-diagonal structure of the factors B±, as well as the expressions for the entries 
of B_, follow from the tri-diagonal structure of the matrix T. The recurrent relation (|8.15| ) 



for the entries of -B+ is equivalent to Bj^B_^ = I + hT. The equations of motion ( pj.l4D are 
now nothimg but the componentwise form of the matrix equation B^T = TB^. ■ 
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The map (|8.14| ) will be denoted dTL. Due to the asymptotics 



^.16|) it is easy to see that 
the equations of motion ( |8.14| ) of the dTL serve as a difference approximation to the Toda 
flow TL The construction assures numerous positive properties of this discretization: 

the map dTL is Poisson with respect to each one of the Poisson brackets (|8.3| ), ( ^.5|) , (|8.7|) , 
it has the same integrals of motion as the flow TL, the Lax representation for the dTL lies 
in its very definition, etc. The only unpleasant property of the equations of motion (|8.14|) , 
when compared with their continuous time counterparts ( ^.1| ) , is the nonlocality. The source 
of nonlocality are the functions (3k- In the open-end case they have explicit expressions in 
terms of finite continued fractions: 



/3fc = 1 + hbk 



1 + hhk-i - 



hbi 



In the periodic case for small h the (3k s may be expressed as analogous infinite (periodic) 
continued fractions. 



8.4 Local equations of motion for dTL 

Fortunately, there exist different ways to bring these equations of motion into a local form 
connected with a localizing change of variables in the sense of Sect. |^. Consider another 
copy of the phase space T. The coordinates in this another copy will be denoted by a^, 
bfc. The localizing change of variables for dTL is the map T(a, b) T{a, b) defined by the 
following formulas: 

Ofe = afc(l + /ibfc), bk = hk + hak-i (8.19) 
The implicit functions theorem assures that this is a local diffeomorphism between the two 



copies of T by small enough values of h. This change of variables appeared in ||K1|| as a Miura 
map in connection with the problem of deformation of integrable systems, but without any 
relation to integrable discretizations. 

Theorem 8.2 The change of variables ( ^.191) conjugates the map dTL with the map 
(a, b) (—> (a, b) described by the following equations of motion: 

afc(l + /ibfc) = afc(l + /ibfc+i), b*. = b*. + /i(afc - afc_i) (8.20) 

Proof. The key point is the following observation: the auxiliary functions (3k aquire in the 
new coordinates local expressions, namely, 

Pk = l + hhk (8.21) 
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To demonstrate this, it is enough to notice that from (|8.19| ) there follows that the quantities 

rrent relation as th 

1 + hhk = 1 + hbk 



1 + hhk satisfy the same recurrent relation as the quantities Pk, namely 



Due to the uniqueness of the solution with the asymptotics 1 + 0{h) the formula ( |8.21| ) is 
proved. The statement of the theorem follows now immediately from ( ^.141) , ( p.21|) . ■ 

Of course, the map ( |8.20| ) is Poisson with respect to pull-backs of the three invariant 
Poisson structures of the Toda lattice. These pull-backs are described by highly nonlocal 
and non-polynomial formulas. However, there exist certain linear combinations of the basic 
Poisson structures in coordinates (a, b) whose pull-backs to the coordinates (a, b) are local. 

Theorem 8.3 a) The pull-back of the bracket 

{•,-}i + M-,-}2 (8.22) 
on T{a,b) under the change of variables (|8.19| ) is the following bracket on T(a, b).• 
{bfc, a^} = -afc(l + /ib^), {afc,bfe+i} = -afc(l + /ibfc+i) (8.23) 

b) The pull-back of the bracket 

{•,-}2 + M-,-}3 (8.24) 
on T{a,b) under the change of variables ( |8.19| ) is the following bracket on T(a, b).- 

{bfc, a^} = -a-kihk + ha.k){l + hhk), {a^, b^+i} = -afc(bfe+i + hak){l + hhk+i) 
{afc, afc+i} = -afcafc+i(l + hhk+i), {b^, b^+i} = -afc(l + hhk){l + hhk+i) 

(8.25) 

c) The brackets ( ^.23|) , ( p.25| ) are compatible. The map (|8.2CI|) is Poisson with respect to both 
of them. 

Proof. To prove the theorem, one has, for example, in the (less laborious) case a) to verify 



the following statement: the formulas ( |8.23| ) imply that the nonvanishing pairwise Poisson 
brackets of the functions ( |8.19| ) are 

{bk, flfc} - afc(l + hbk), {ak, bk+i} = -afc(l + hbk+i) 

{a/c5«fc+i} = —hakak+i, {bk,bk+i} = —hak 

This verification consists of straightforward calculations. In what follows we do not repeat 
analogous arguments in the similar situations. I 
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The map ( p.2(]| ) was first found in ||HTI|| , along witli tlie Lax representation. In [ po| it was 
stressed tliat tfiis map is notliing otlier tlian tlie so-called qd algorithm well known in the 
numerical analysis. Its Poisson structure and its place in the continuous time Toda hierarchy 
were not discussed in [ HTl ]. The previous theorem provides a bi-Hamiltonian structure of 
the qd algorithm. This result in a slightly different form was found in 



Theorem 8.4 |KT]. The pull-hack of the flow TL under the change of variables 
described by the following differential equations: 



1|) is 



kk = afc(b 



k+l 



(afc - afe_i)(l + hhk 



(8.26) 



Proof. To determine the pull-back of the flow TL, we can use the Hamiltonian formalism. 



An opportunity to apply it is given by the Theorem |8.3| . We shall use the statement a) only. 
Consider the function /i~^Hi(a, 6) = Y,k=i^k- It is a Casimir of the bracket {t}!, and 
generates exactly the flow TL in the bracket h{-, ■}2. Hence it generates the flow TL also in 



the bracket ( 8.22 ). The pull-back of this Hamilton function is equal to J2k=ii^k + hak-i) ■ 
It remains only to calculate the flow generated by this function in the Poisson brackets ( |8.23D . 
This results in the equations of motion ( p.26|) . ■ 



9 Second flow of the Toda hierarchy 

Now we want to demonstrate that our method for finding integrable discretizations and local 
equations of motion for them works not only for the flow TL, but equally well for the higher 
flows of the Toda hierarchy. We consider here the second flow (called hereafter TL2). 



9.1 Equations of motion and tri— Hamiltonian structure 

This is the flow on T governed by the differential equations 

Ofc = afc(&fc+i - + ctfe+i - ctfc-i), bk = ak{bk+i + bk) - ak-i{bk + bk-i) (9.1) 



This flow is Hamiltonian with respect to all three brackets 
sponding Hamilton functions are: 



(^.7|). The corre- 



H3(a, b) = -J2bl + J2 Uak + ak-i) = -{tiT^a, b, A))^ 



(9.2) 



k=l k=l 

for the bracket {■, B.2{a, b) for the bracket {■, ■}2, and Hi(a, h) for the bracket {•, - ja. 
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9.2 Lax representation 

Naturally, the Lax matrix for the flow TL2 is the same as for the flow TL. The difference 
lies in the auxiliary matrices B± taking part in the Lax representation (since in this section 
we are dealing only with the flow TL2, using the same notations B± as in Sect. |^ will not 
lead to confusions; the same holds also for some other notations in this section). 
The equations of motion ( p.l| ) are equivalent to the Lax equations in g 



so that 



f = [T, B+] = -[T, 5_] with 5± = 7r±(T2) (9.3) 



N N N 

B+{a,b,\) = Y.^hl + ak + ak-i)Ek,k + \Y.^hk+i + hk)Ek+i,k + >^^Y.^k+2,k (9.4) 

fc=l k=l k=l 

N N 

5_ (a, 6, A) = A"^ ^ {hk+i + hk)akEk^k+i + A"^ ^ ak+iakEk^k+2 (9.5) 

k=l k=l 

9.3 Discretization 

In order to obtain an integrable discretization of the flow TL2 we can apply the recipe of 
Sect. ^ with F{T) = / + /iT^, i.e. consider the map described by the discrete time Lax 
equation 

f = B-^TB^ = B_TB_^ with B± = n±(J + hT"^) (9.6) 

Theorem 9.1 The discrete time Lax equation ( |9.(j| ) is equivalent to the map {a,b) t— > (a, 6) 
described by the following equations: 

Pk+i r , , ,fo'kSk ak-iSk-i\ „N 
flfc = flfc —5—, bk = bk + h\ — (9.7) 

Pfc V Pfc Pfc-i / 

where the auxiliary functions 6k = 6k{a, b) = 0(1) and [3k = Pk{o,, b) = 1 + 0{h) are uniquely 
defined for h small enough by the recurrent relations 

Sk = bk^. + bk- ^^^^^ (9.8) 
Pk-i 

n 1 , L/t2 , , ^ ^^O-fc-lf^fc-l h'^0,k~ia'k~2 /^i 

I3k = l + h{bk + ak + ak-i) (9.9) 

Pfc-l Pk-2 

and have the asymptotics 

5k = bk+i + bk + 0{h) (9.10) 

(5k = l + h{bl + ak + ak-i) + 0{h'') (9.11) 
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Remark. The matrices B± from this theorem have the following expressions: 

N N N 

B+ia, b,X) = J2 PkEk,k + h\J2 hEk+i,k + ^ Ek+2,k (9.12) 

k=l k=l k=l 

B.{a, b,X)=I + hX-' ^E,,,^i + hX-' E ^^E,,,^, (9.13) 

k=i k=i 

Proof. The tri-diagonal structure of the matrix T assures that the factors B± have the 
structure as in (|9.12|) , (|9.13|) . The expressions for the entries of B^ and the recurrent 



relations for the entries of -B+ follow easily from the equality B+B^ = I + hT"^. After that 
the equations of motion follow from the equation B^T = TB^. ■ 

9.4 Local equations of motion for dTL2 



For the map (|9.7|) , called hereafter dTL2, there exists a localizing change of variables, however 
different from the one used for the map dTL. It is given by the formulas 

ak = SLk{l + hak-i){l + hhl), bk = hk{l + hak-^i) + hhk-iak-i (9.14) 

(Actually, discretizations of all higher flows of the Toda hierarchy should possess their own 
localizing changes of variables). 

The name "localizing change of variables" is justified by the following theorems. 

Theorem 9.2 The change of variables ( p. 141 ) conjugates the map dTL2 with the map 
(a, b) I— > (a, b) governed by the following local equations of motion: 

afe(l + /i5fe-i)(l + hhl) = afe(l + hak+i){l + hhl^^) 

(9.15) 

bfc - bfe = hak(hk + hk+i) - ha.k-i(hk-i + hk) 

Proof. This time the key point of the proof is obtaining the following local expressions for 
the coefficients of the factor 11+ (/ + hT"^): 

13k = (l + /iafe)(l + /iafe_i)(l + /ib2) (9.16) 
4 = (l + /ia,)(bfc + bfc+i) (9.17) 



Indeed, the equations of motion ( |9.15|) follow directly from ( |9.7|) and the latter formulas. 



To prove the latter formulas, define the quantities /3k, Sk by the equations ( |9.16| ), ( |9.17| ). 
A straightforward calculation convinces that then the recurrent relations (|9.9|), (|9.8|) are 



22 



satisfied. The uniqueness of solution to tliese equations proves tlie desired expressions for 

Tliere still exists a linear combination of invariant Poisson brackets of the Toda hierarchy, 
whose pull-back with respect to the above change of variables is local. 



Theorem 9.3 The pull-back of the bracket 

{■,-}i + M-,-}3 



(9.18) 



on T{a, h) under the change of variables ( |9.14| ) is the following bracket on T(a, b) 

{bfc,afc} = -afc(l + /iafc)(l + /ib^) 

{afc,bfe+i} = -afc(l + /iafc)(l + /ib^_^i) 
The map ( p.l5| ) is Poisson with respect to ( |9.19| ). 
Proof - by straightforward verification. ■ 



(9.19) 



Theorem 9.4 The pull-back of the flow TL under the map (|9.14|) is described by the fol- 
lowing differential equations: 



kk = afe(bfc+i - bfc)(l + hak) 



(afe - afc_i)(l + hh 



(9.20) 



Proof. The equations of motion we are looking for describe the flow TL as a Hamiltonian 
system in the Poisson bracket (|9.18|) . The corresponding Hamilton function is, obviously, 
given by 

(/ + /iT^) VH(T) = T 

Hence we have H(T) = (2/i)-i(logdet(/ + /iT^))^. From the expressions for the factors of 
the generalized LU factorization of the matrix / + hT'^ given in the formulas (|9.12|) , (|9.13| ) 
we conclude that H(T) = (2/i)^^ X^^Li log(/5fc)- Taking into account the expressions (|9.16| ), 
we flnd flnally: 



1 



N 



N 



k=l k=l 



(9.21) 



The Hamiltonian equations of motion generated by this Hamilton function and the Poisson 
brackets ( p.l9| ), coincide with ( |9.20|) . ■ 

The differential equations ( p.20| ) are a particular case of a two-parameter deformation of 



TL found in |Kl|]. 
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10 Volterra lattice 



10.1 Equations of motion and bi— Hamiltonian structure 

The second flow ( |9.1|) of tlie Toda liierarcliy, unlike the first one, allows an important reduc- 
tion bk = 0, and it results in the famous Volterra lattice (hereafter VL), known also under 
the names of Lotka- Volterra system, discrete KdV equation, Langmuir lattice, Kac-van 
Moerbecke lattice, etc. Q, 



o-k — a-ky^k+i — ctfc-ij 
The phase space of VL in the case of periodic boundary conditions is 

V = M^(ai, ...,aN) 



(10.1) 



(10.2) 



Clearly, it is the subspace 6^ = of T(a, b). Unfortunately, neither of three invariant Poison 
brackets {■, {■, ■}2, and {■, - ja of the Toda hierarchy can be properly restricted to the set 
V. However, it is easy to see that the quadratic bracket {•, ■}2 allows a Dirac reduction to 
this set, the reduced bracket being defined by the relations 

{ak,ak+i}2 = -atak+i (10.3) 
The system VL is Hamiltonian with respect to this bracket with the Hamilton function 



Hi (a) 



N 
k=l 



There exists one more local Poisson bracket on V invariant with respect to VL and compatible 
with ([HOI) 



-O.kO.k+lO.k+2 



(10.4) 



FTTp. It is given by the relations 

{ofc, Ofc+ija = — flfcO.A;+i(afc + Ofc+i), {flfc, 0^+2)3 
the corresponding Hamilton function being equal to 

1 ^ 

Ho(a) = -J^logK) 
^ k=i 



10.2 Lax representation 

The Lax representation of the previous section survives in the process of restriction to V. 
So, we obtain the Lax representation of the VL with the matrices from g 

f = [T,B^] = -[T,B_] 



KM| 



:io.5) 
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where 

N N 



T{a,X) = X-^J2('kEk,k+i + Xj2Ek+i,k (10.6) 

fc=l k=l 

N N 

B+ia,\) = 7r+(T2) = ^(a^ + afc„i)Efc,,, + A' J] Efc+2,fe (10.7) 

k=l k=l 

N 

B.{a,X) = 7r_(T2) = A'^ ^ afc+iafcEfc,fc+2 (10.8) 

k=l 

These equations allow an r-matrix interpretation in case of the quadratic bracket {■, ■}2, see 
||S4|| . The Volterra hierarchy consists of the flows allowing Lax representations of the form 
( |10.5| ) with B± = 7r±(/(T^)), where f : g ^ g is Ad*-covariant. 

10.3 Discretization 

The discretization of the second Toda flow ( p. 71) also may be restricted to the set V. It is 
easy to see from ( p.8| ), ( |9.7| ) that in this situation we have with necessity 6k = 0, and the 
recurrent relations (|9.9| ) turn into 

A = 1 + h{ak + a,_i) - (10.9) 

Hk-2 



Theorem 10.1 []S^ (see also [[K2|| ). The discrete time Lax equation 



T = B-^TB+ = B TBZ^ with B± = n±(/ + hT^) (10.10) 
is equivalent to the following map a ^-^ a: 

ak = ak^^ (10.11) 

where the quantities jSk = /3fc(a) = 1 + 0{h) are uniquely defined by the recurrent relations 

(3k - ha, = ^ ^^-] = I + ^ (10.12) 
Pk-i — n-o-k-i Pk-i — n-o-k-i 

and have the asymptotics 

pk = l + h{ak + ak-i)+0{h^) (10.13) 
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Remark. The matrices B± have the following expressions: 

N N 

B+ia,X) = U+il + hT') = Y.PkEk,k + h\^Y.Ek+2,k (10.14) 

fc=l k=l 
N 

B.{a,X) = 1141 + hT') = I + h\-'Y.^^E,,,^2 (10.15) 

k=l Pk 



Proof. Referring to Theorem 9.1, we need only to prove that the auxiliary quantities 



f3k defined by the recurrent relations ( |10.9| ), may be alternatively characterized as solutions 



to (|10.12|) . To do this, consider the matrix equation B+T = TB^. In coordinates it is 
equivalent to the set of two scalar equations: ( p.4.17| ) and 



Pk-l 

The last equation means that the following quantity does not depend on k: 

{(3k - hak){(3k-i - hak^i) 



c = const (10.16) 



Pk-i 

and we need only to prove that c = 1. But applying ( |10.16| ) twice, we have: 



n u , hcak-i , hak-i{(3k-2 - hak-2) , , h^ak-iak-2 
Pk - httk = c + = c H = c + httk-i ^ 

Pfc-l — iia-k-l Pk-2 Pk-2 

Comparing this with (|10.9| ), we conclude that c = 1. ■ 



The map ( p.0.11| ), ( |10.12| ) will be denoted dVL. It is Poisson with respect to the brackets 
( |10.3| ), (|10.4|) , but nonlocal, as opposed to the flow VL. In the open-end case we have 
terminating continued fractions for the nonlocal quantities [3k (in the periodic case for small 
h the Pk may be expressed as analogous infinite periodic continued fractions). From ( |10.9| ) 
we derive, for example, for l3k with odd indices: 

/?2fc+i = 1 + h{a2k+i + a2k) - -7—n ^ '^2fcQ2fc-i (10.17) 

l + h{a2k-i + a2k-2)- h'^a^ai 

1 + hai 



(with obvious changes for j32k)- Analogously, from ( |10.12| ) we derive alternative continued 
fractions: 

I3k~hak = l + ^^^T^, (10-18) 

1 + • . . _^ ha2 



1 + hai 
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One says that the continued fractions (|10.17|) are obtained from ( p.0.18| ) with the help of 
compression. 



10.4 Local equations of motion for dVL 



The locahzing change of variables V(a) i-^ V(a) for dVL follows from ( |9.14| ) upon setting 
bfe = 0: 

ak = akil + hak-i) (10.19) 



Again, this change of variables was found in ||K1|| in the context of integrable deformations, 
but with no relation to integrable disacretizations. 



As it follows from the proof of Theorem |9]^, in the variables the auxiliary quantities 
Pk become local: 

/3fc = (l + /iafe)(l + /iafe_i) (10.20) 
and we obtain the following result. 

Theorem 10.2 The change of variables ( p.0.19| ) conjugates the map dVL with the map a i— 
a governed by the following local equations of motion: 

afc(l + /lafc-i) = afc(l + /lafc+i) (10.21) 



The only known local invariant Poisson bracket ( |9.19| ) for the local form of dTL2 does 



not admit a proper restriction or a reduction to the subset V(a) of T(a, b) characterized 
by bfc = 0. Nevertheless, there appears to exist a local bracket on V(a) invariant under the 
local form of dVL. 

Theorem 10.3 The pull-back of the bracket 

{■,-}2 + M-,-}3 

on V(a) under the change of variables ( |10.19 ) is the following bracket on V(a).- 



{ak, afc+i} = -afcafc+i(l + hak){l + hak+i) (10.22) 



This bracket is invariant under the map (|10.21 



Proof - by direct calculation. ■ 

The local equations of motion (|10.21D together with Lax representation were found by 



different methods in [irHOII, [|FN||, |E31. The invariant Poisson structure and the relation to 



factorization problem were pointed out only in the latter of these references. 
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Theorem 10.4 ||K1|| . The pull-back of the flow VL under the map ( p.0.19| ) is described by 
the following equations of motion: 



kk = ak{l + hak){ak+i - ak-i) 



[10.23) 



Proof. Since the flow VL has a Hamihon function (2/i)~^ J2k=i log('^/c) the Poisson bracket 
h{-, - js, and this function is a Casimir of the bracket {•, •}2, we conclude that this function 
generates the flow VL also in the bracket {■, ■}2 + h{-, ■}^. That means that the pull-back 
of the flow VL is a Hamiltonian flow in the bracket ( |10.22| ) with the Hamilton function 



N 



(2/i)-^^log(afc(l + /iafc 



k=l 

Calculating the corresponding equations of motion, we arrive at ( 10.23 ). ■ 
The system ( |10.23 ) is known under the name of the modified Volterra lattice. 
Remark. An additional map 

afc ^ Ofc = ——T — 
1 + hak 

delivers an alternative version of the localizing change of variables 

_ Ofc 

(1 - hak){l - httk-i) 

The corresponding local version of the dVL reads: 



;i0.24) 



(1 - hak){l - hak- 



1 - hak){l - hak+i) 



;i0.25) 



This version was introduced in [[HT|| under the name "discrete Lotka- Volterra equation of 
type II" (while the type I was assigned to ( |10.21|) ). The modifled Volterra lattice in the 
variables takes the following form: 



1 - /iflfeH 



1 - hak-i 



;i0.26) 



A remarkable feature of the variables (not mentioned in [[HT|| ) is the formal coincidence 
of the invariant Poisson bracket ( |10.22[ ) with the quadratic invariant Poisson bracket of the 
continuous time flow VL. Indeed, in the variables the bracket (|10.22| ) takes the form 

{ak, ak+i} = -akak+i (10.27) 
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11 Second flow of the Volterra hierarchy 



Now we apply the general procedure of integrable discretization to the second flow of the 
Volterra hierarchy. We use the notations /3k, B± etc. for objects analogous to those of the 
previous section without danger of confusion. 

11.1 Equations of motion and bi— Hamiltonian structure 

The second flow of the Volterra hierarchy (hereafter VL2) is described by the following 
differential equations on V: 

ak = ak(ak+i{ak+2 + ak+i + cik) - afc-i(afc + ctfc-i + ak-2)) (H-l) 



This flow is Hamiltonian with respect to both Poisson brackets ( 10.3 ), ( 10.4 ). The corre- 
sponding Hamilton functions are: 

N N 

k=i k=i 

for the quadratic bracket {■, ■}2, and Hi(a) for the cubic bracket {■, - js. 
11.2 Lax representation 

The Lax representation for the flow VL2 is of the type (|3.1]) with /(T) = T^. 
Theorem 11.1 The flow (|1 1 . 1|) admits the following Lax representation in g: 

t=[T,B+] = -[T,BJ\ (11.2) 

with the matrices 

N 

B+{a,X) = TT+{T^) = ^ (cLk+iak + {ak + ak^i^ + ak^iak-2)Ekk 

k=l 

N N 

+X ^^{ak+2 + ttk+i + + ak-i)Ek+2,k + X Ek+i,k 

k=l k=l 

N 

B_{a,X) = 7T_{T'^) = X'"^ '^{ak+2 + cik+i + cik + ak-i)ak+iakEk,k+2 

k=l 

N 

+ X ^ Z 0,k+30'k+2C'k+lO'kEk,k+4: 
k=l 
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11.3 Discretization 

Applying the recipe of Sect. ^ with F{T) = I + hT'^, we take as a discretization of the flow 
VL2 the map described by the discrete time Lax equation 

f = B-^TB+ = B TBZ^ with S± = n± (j + /iT^) (11.3) 



Theorem 11.2 The discrete time Lax equations ( p.1.3 ) are equivalent to the map a ^ a 

described by the following equations: 

Sfc = Ofc %^ (11.4) 

Pk 

where the functions Pk = Pk{ci) = 1 + 0{h) are uniquely defined for h small enough simulta- 
neously with the functions 6k = 0(1) by the system of recurrent relations 

Pk - h{6k - ak+2)ak = uf/^'^ ^ (11-5) 

Pfc-i — iiyOk-i — Ofc+ija,fc-i 

^fc = afc+2 + ctfc+i + ctfc + flfc-i n (11-6) 

Pfc-2 

The auxiliary functions Pk have the asymptotics 

Pk = l + h(ak+iak + {ak + ak-i)"^ + ak-iak-2) + 0(/i^) (11.7) 
Remark. The matrices B± have the following expressions: 

N N N 

B+{a,\) = ^PkEk,k + h\'^^5kEk+2,k + h\^^Ek+4^k 

k=l k=l k=l 

N x N 

B4a,X) = I + hX-^ E ^^^Ek,k-,2 + hX-^ E Ek,k^, 

k=l k=l 

Proof. The scheme of the proof is standard. First of all, the general structure of the 
factors B± is clear from the bi-diagonal structure of the matrix T. The expressions for 
the entries of B_ and the recurrent relations for the entries of -B+ follow from the equality 



B^B_ = I + hT^. However, in this way we come to the recurrent relation (|11.6|) for 6k and 
the following recurrent relation for Pk'. 

n . , if , / , \2 , \ h'^ak-iak-25l_2 h'^ak-iak-2ak-^cik-A „x 
Pk = l + h[ak+iak + [ak + ak-i) +afc_iafc_2 (11.8) 

^ ^ Pfc-2 Pfc-4 
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different from (|11.5|) . In order to "decompress" this recurrent relation for and to bring it 



into the form ( |11.5D , we proceed as follows. The matrix equation B^T = TB^ is equivalent 



to the set of three scalar ones: the equation of motion (|11.4| ) and 



Pk+i - hSkttk+i = iPk-i - hSk^ia^-i) (11.9) 

Pk-l 

5k+i - ak+3 = 5k- CLk-i (11.10) 

Pfc-i 

From the last two equations it is easy to derive: 

Pk+i — h{5k+i — ak+sjttk+i = -w~ (Pk-i — h{5k-i — ak+i)ak-i) 

Pk-i ^ ^ 

which means that the following quantity does not depend on k: 

(Pk - h{5k - ak+2)ak) (Pk-i - h{5k-i - ak+i)ak-i) 

^ = c = const (11-11) 

Pfc-i 

It remains to prove that this constant is equal to 1. To this end we apply ( |11.11D twice (as 



in the proof of Theorem |10.1| ) to derive 

h^{5k-l — 0'k+l){Sk-2 — (^k)(^k-lO'k-2 



i3k = c + h{5k - ak+2)ak + K^k^i - ak+i)ak-i 



k-2 



Transforming the second and the third terms on right-hand side with the help of ( |11.6| ) 



and the fourth one with the help of ( p.l.lO| ), we see that the 0(/i)-terms of the last formula 



exactly coincide with the 0(/i)-terms on the right-hand side of ( |11.8D . Equating the free 
terms, we get c = 1. ■ 

11.4 Local equations of motion for dVL2 

The localizing change of variables for the map dVL2 is given by the formulas 

_ (1 + hal^^) 

(1 - /iafcafe_i)(l - /iafc_iafc_2) 

Theorem 11.3 The change of variables ( |11.12|) conjugates the map dVL2 with the map 
a I— > a described by the following local equations of motion: 

{l + hal_,) _^ (l + /^ag^i) 

(1 - /iafcafe_i)(l - hak-iak-2) (1 - hak+2ak+i){l - hak+iak) 
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Proof. The statement will follow immediately, if we prove the following expressions for the 
auxiliary quantities Pk- 

= 71 7 -\n 7^ ^2n . 7 11-14 



To do this, it is sufficient to verify that the recurrent relations (|11.5|) , (|11.6|) are satisfied by 
the quantities ( |11.14| ) and 

^ + a^^i ^ a^ + afc_i 

(1 - hak+2ak+i){l - hak+iak) (1 - hak+iak){l - hakak-i) 

Such a verification is a matter of straightforward, though somewhat tedious algebra. We 
omit it, giving only an important intermediate formula: 

- h{6k - a,^2)ak = jz ^'^ \^^'-]! r (11.16) 

(1 - /iafcafc_i)(l - hak-iak-2) 

(it implies immediately (|11.5|) ). ■ 

Unlike the situation with the previously encountered localizing changes of variables, we 
did not find a linear combination of invariant Poisson brackets of the Volterra hierarchy, 
which would be pulled back into a local bracket on V(a) (presumably, such a linear com- 
bination does not exist). Therefore, to find a pull back of the fiow VL under the change 
of variables ( |11.12| ), we cannot use the Hamiltonian formalism and are forced to turn to a 



direct analysis of equations of motion. 



Theorem 11.4 The pull-back of the flow VL under the change of variables (|11.12|) is de- 
scribed by the following equations of motion: 

' ''^ ^ ^Hl-/iafc+iafc l-hakak-ij (H-l"^) 



Proof. It is a matter of straightforward calculations to verify that the equations ( |11.17| ) are 
sent into ( |10.1|) by the change of variables ( |11.12| ). ■ 

11.5 Local discretization of the KdV 

Probably the most interesting feature of the fiow VL2 is that a certain linear combination of 
this fiow with the original VL gives a direct spatial discretization of the famous KdV | |AL 1|| . 
Indeed, the equations of motion of the linear combination a ■ VL + VL2 read: 



(afc+i(a + afc+2 + Ofc+i + Ofc) - afc-i(a + + a^-i + 0^-2)) (11.18) 



32 



Setting a = —6 and 



afc = 1 + e'^pk 



we get the equations of motion 

Pk = (l+e'^pk) {pk+2 - 2pfc+i + 2pk-i -pk-2+£^Pk+i {Pk+2 +Pk+i +Pk) -e^Pk-i{Pk +Pk~i+Pk-2)) 

Assuming that Pk{t) p(t, ke) with a smooth x) and rescahng the time t ^ t/ (2£^), we 
see that the previous equation approximates 

Pt = Pxxx + 'iPxP 

which is the KdV. Therefore a local discretization of the linear combinations of the flows VL 
and VL2 will result (for a = —6) in a local spatio-temporal discretization of the KdV (see 
r^] for a highly nonlocal integrable discretization). 



Theorem 11.5 A localizing change of variables for the discretization of the above flow cor- 
responding to F{T) = I + h{aT'^ + T^) is given by the formulas 

(1 + haak-i + /la^.J (1119) 
(1 - /iafeafc_i)(l - hak-iak-2) 

This change of variables conjugates the above discretization with the map a 1-^ a described 
by the local equations of motion 

(1 + haak-i + hal_^) _ (1 + haak+i + hal^^) 

(1 - /iafcafc_i)(l - hak~iak-2) (1 - hB.k+2ak+i){l - hak+i^k) 



Proof is completely analogous to that of Theorem 11.3, and will not be repeated here. 



12 Modified Volterra lattice 

12.1 Equations of motion and Hamiltonian structure 

As we have seen many times in the preceding exposition, the localizing changes of variables 
are interesting not only in their connection with the problem of integrable discretization, 
but already on the level of continuous time systems. Namely, they may be viewed as Miura 
transformations and used to find the so called modified equations of motion. These modified 
systems are often interesting in their own rights - and one can also wish to discretize them 
as well. We consider in the present section one example of modified systems, namely, the 
modified Volterra lattice (MVL). 
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To this end we use the change of variables (|10.19|) in shghtly different notations. Namely, 
we define a Miura change of variables A4 : A4V{q) ^— V(a) by the formula 

ak = qk{l + aqk-i) (12.1) 

We denote the parameter a instead of h {a is not related to the time step of discretizations 
and is not supposed to be small). The MVL is the pull-back of the VL (|10.1| ) under the 
change of variables ( |12.1| ). Its equations of motion were already determined in Sect. |T^ to 
be 

gfc = gfc(i + agfe)(gfc+i - gfc-i) (12.2) 

As usual, this system may be considered under periodic or open-end boundary conditions, 
the phase space in case of the periodic ones being 

A<V = M^(gi,...,g^) 

As pointed out in Sect. this system is Hamiltonian with respect to the following Poisson 
bracket on A^V: 

Uk, Qk+i} = -QkQk+iil + aqk){l + aqk+i) (12.3) 
and with one of the Hamilton functions 

N N 

Ho(g) = 5:iog(gfc) or Ho(g) = ^ log(l + ag^) (12.4) 

fc=l k=l 



(their difference is a Casimir of the bracket (|12.3| )). The bracket (|12.3|) is the pull-back 



under the Miura change of variables (|12.1|) of the invariant Poisson bracket {■, ■}2 + «{-, - js 



of the VL (see ( |10.3| ), ( |10.4| ) for the relevant definitions) 



12.2 Discretization 



We define the discretization dMVL of MVL ( p.2.2 ) as a pull-back of the dVL under the 
Miura change of variables (|12.1|) . 



Theorem 12.1 The equations of motion of the map dMVL read: 

qk = qk (12.5) 

Ik 

where the quantities •jk = lk{q) = 1 + 0{h) are uniquely defined by h small enough by the 
system of recurrent relations 

, Ik-i + haqkqk-i ^ , ^ ,^ , ^ ^7^^^ 

lk-hqk = T -ik - hqk + aqk = [I + aqk) r (12.6) 

7fc_i - hqk-i 7yfc-i - hqk-i 

and have the asymptotics 

7fc = 1 + h{qk + qk-i + aqkqk-i) + 0{h^) (12.7) 
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Proof. From the equations of motion (|12.5| ) and the definition (|12.6|) we derive: 

7fc_i (7fc+i - hqk+i] 



l + aqk = — (7fc + haqk+iqk) ( H 

Ik \ Ik-hqk^ 



aqk) 



Ik {ik-i - hqu-ij 

Hence the variables au defined as in ( [L2.1| ) satisfy the equations of motion (|10.11| ) with 

Ik 



12.8) 



(3k 



lk-2 



{ik-i - hqk-i)(cik-2 - hqk-2) 



12.9) 



It remains to prove that these quantities jSk satisfy the recurrent relations ( [10.121 ). But this 
follows immediately from the formula 

Ik-i 



Pk - hak 



lk-2 



{lk-2 - hqk-2) 



which is an easy consequence of (|12.9|) , ( |12.1| ), and (|12.6|) . ■ 

Clearly, by definition the map dMVL is Poisson with respect to the bracket ( 12. 3|) , but 
the above theorem renders dMVL nonlocal. In particular, in the open-end case we have 
terminating continued fractions 

hqk-iil + aqk) 



lk-hqk = l + 



hqk- 



aqk-i) 



;i2.io) 



1 + 



+ 



hq2{l + aq-s) 



1 + hqiil + aq2) 
In the periodic case these continued fractions are also periodic. 

12.3 Local equations of motion for dMVL 

It is possible to find a localizing change of variables A^V(q) ^ M.V{q) for the map dMVL: 

* = 3i(l±^ (12.11) 

1 - /iaqfcqfc-1 

As usual, this map is a local diffeomorphism for h small enough. 

Theorem 12.2 The change of variables ( |12.11|) conjugates the map dMVL with the map 
q 1-^ q described by the following local equations of motion: 



qfc 



1 + h<\k-i) 



qfc 



or, equivalently, 



1 + «qfc 1 + aqfc 

1 + aqfe 1 + aqfc 



(1 + hc^k+i) 



1 - /laq^qfc-i 1 - /laqfe+iq^ 



;i2.12) 
;i2.13) 
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Proof. The first step is to prove the following local expressions for the quantities 7^ in the 
variables q^: 

7. = (i + ^q^)(i + ^q^-i) (12.14) 

1 - haqk(lk-i 



Indeed, from ( [I2.5| ), (|12.11|) and (|12.14|) we derive immediately the equations of motion in 
the form 

qA;(l + hqk-i) _ qfc(l + hqk+i) 
1 - haqk^k-i 1 - haqk+iqk 

which is equivalent to either of (|12.12|) and (|12.13|) . To prove ( p.2.14| ), it is as usual enough 
to verify that the quantities defined by this formula satisfy the recurrent relations (|12.6|) . 
This verification is the matter of a simple algebra. ■ 

The equations of motion ( |12.12| ) were introduced in | pd'l 



In general, the only expression for a Poisson bracket invariant with respect to the map 
( p.2.12| ) is non-local and non-polynomial (it may be characterized as the pull-back of the 
bracket (|12.3 ) under the localizing change of variables ( 12.11|) ). However, a direct analysis 



of equations of motion allows to prove the following statement. 



Theorem 12.3 The pull-back of the flow ( |12.2|) under the change of variables ( |12.11|) 
described by the following equations of motion: 



IS 



qfc = q,(l + aq,)(l + /iq,) ( -—^ ^-^^ ) (12-15) 

V 1 - haqk+iCLk 1 - /iaq^qfe-i / 

This system may be called the second modification of the Volterra lattice. Notice its remark- 
able symmetry with respect to the interchange a ^ h. Also the expressions for in terms 
of c[k enjoy this symmetry. Indeed, the composition of two changes of variables (|12.1| ) and 
( I12T1I ) resilts in 

qfc(l + /iqfc_i)(l +aqfc_i) 

= n h vn h V (12.16) 

(1 - toqfcqA:-i)(l - haqk-i(lk-2) 



However, in the discrete equations of motion ( |12.12| ) this symmetry goes lost. 



The formula ( |12.16| ) allows also to translate the Miura map A4 : A4V{q) ^ V(a) into 
the language of localizing variables. Namely, if we define the map M : AiV{q) ^ V(a) by 
the formula 

a. = 3i(l±i^3tl) (12.17) 

1 - haqkq,k-i 

then an easy calculation shows the commutativity of the following diagram: 
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M 

MV{q) V(a) 



dmri) 



MV{q) 



V{a) 



M 



Remark. Sometimes another form of MVL is more convenient: 

Ck = {cl-e^){ck+i-Ck-i) (12.18) 

It is related to (|T23) with a = (26) by means of a hnear change of variables: 

qk = Ck-e (12.19) 

and time rescaling t ^ t/a. Analogous linear change of localizing variables = — e 
accompanied by the rescaling of the time step h \—>- h/a allows to derive from ( |12.12| ) the 
following local integrable discretization of the flow ( [12.181 ): 



Ck-e 



Ck + e 

12.4 Particular case a 



1 + 2/i£:Cfc_i) 



2hec 



^12.20) 



oo 



There exists a particular case of MVL, for which a more detailed information is available. It 
corresponds to the a ^ oo limit of the previous constructions, which has to be accompanied 
by the time rescaling t ^ t/a. (Alternatively, one could consider the e limit of the 
system ( |12.18|) ). Let us list the corresponding results. The equations of motion read: 



Qk = ql{qk+i -qk- 



^12.21) 



and the Miura map 7V1 : M.V{q) ^ V(a) relating this system to the VL, takes the form 

ak = qkqk-i (12.22) 
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An invariant Poisson bracket ( |12.3|) in the limit a — > cxd turns (being suitably rescaled) into 

{qk, gfc+i}3 = -qM+i (12.23) 

(see below for the reason for assuming the index "3" to this bracket). The Hamilton function 
of the flow (|12.21|) in this bracket is equal to 



TV 

Ho(g) = ^log(gfc) 

fc=i 



The localizing change of variables ( |12.11| ) after rescaling h \^ h/a and going to the limit 
a — > cxD becomes 

= (12.24) 

and the discrete time equations of motion ( |12.13| ) become 



;i2.25) 



1 - /iQfcqfc-i 1 - /^qfcqfe+i 
or, equivalent ly, 

q/c - q/c = hc^k(ik{(ik+i - qfc-i) (12.26) 

Finally, the pull-back of the system ( p.2.21| ) under the localizing change of variables ( |12.24| ) 
reads: 

n,, = n? I 

1 - /iqfc+iqfc 1 - /iq^q/c-i 



^^ = ^k[-r^ . (12.27) 



Now notice that the bracket ( p.2.23| ) is the pull-back with respect to the Miura map M. 
( p.2.22| ) of the cubic invariant bracket {■, -js of the Volterra hierarchy. It turns out that in 
this particular case the pull-back of the quadratic bracket {•, •}2 is given by more elegant 
formulas than in the general case. Namely, the corresponding bracket on }AV[(i) is quadratic 
albeit still nonlocal: 

{qk,qj}2 = T^kjQkqj (12.28) 

with the coefficients 



k-j=0 
^fcj ~i^l k — j > and odd or k — j < and even 

1 k — j > and even or /c — j < and odd 



;i2.29) 



The corresponding Hamilton function of the flow ( |12.21|) is given by 



H(?) = J2 qkqk+1 

k=l 
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So, in this case we know a simple formula for an additional invariant Poisson bracket. This 
allows to find a nice formula for an invariant Poisson structure of the map ( |12.26D and of 
the flow (|12.27|) - the result lacking in the general case. 



Theorem 12.4 The pull-back of the bracket 

{■,-}2 + h{;-h 

on AiV{q) under the map (|12.24|) is the following bracket on A4V{q): 



;i2.30) 



^12.31) 



with the coefficients ( |12.29|) . The map (|12.26|) and the flow ( |12.27] ) are Poisson with respect 
to this bracket. 



Proof. A direct calculation shows that the brackets ( |12.31|) for the variables are sent by 
the map ( |12.24|) into the brackets (|12.3CI|) for the variables qk- I 



We close this section by noticing that the translation of the Miura map Ai ( |12.22| ) into 
the localizing variables is the map M : AiV{q) ^ V(a) given by the formula 



1 - hc[k<ik-i 

For the localizing variables of VL introduced by ( p.0.24] ), we have: 

= q^q/c-i 



;i2.32) 



;i2.33) 



13 Bogoyavlensky lattices 

13.1 Equations of motion and Hamiltonian structure 

There are three basic families of integrable lattice systems carrying the name of Bogoyavlen- 



sky 



(although some of them were found earlier in 



These systems are 



enumerated by integer parameters m,p > 1 (p > 1 for the third one) and read: 



dfc — ttk 










i=i 




( V 


p 




n ^k+j 


- n 









;i3.i) 



:i3.2) 
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flfc = n - n (13.3) 

We shall call these systems BLl(m), BL2(p), and BL3(p), respectively. 

The lattices BLl(m) and BL2(p) serve as generalizations of the Volterra lattice, which 
arises from them by m = 1 and p = 1, respectively. The lattice BL3(p) after the change of 
variables ak ^ qk = ^ and t ^ —t turns into 

/p-i p-i \ 

= n ^k+j - n (Ik- J (13-4) 

We call the latter system modified BL2(p). It serves as a generalization of the modified 
Volterra lattice ( |12.21|) , which is the p = 2 particular case of ( |13.4D . 

These systems may be considered on an infinite lattice (all the subscripts belong to Z), 
and admit also periodic finite-dimensional reductions (all indices belong to Z/NZ, where N 
is the number of particles). The lattices BLl and BL2 admit also finite-dimensional versions 
with the open-end boundary conditions. The phase space of the periodic BL's is 

i3 = K^(ai,...,ajv) (13.5) 



The Hamiltonian structure of the BLl was determined in Kl and later in [B|; for the BL2 



and BL3 in the infinite setting this was done in ||ZTOF|| ; in the open-end and the periodic 
setting (where some subtleties come out) the Hamiltonian structures were determined in 
S4| . We reproduce here the corresponding result for the periodic boundary conditions. 

The invariant quadratic Poisson brackets for Bogoyavlensky lattices are given by the 
formula (we set p = 1 for BLl(m), m = 1 for BL2(p), and m = —1 for BL3(p)): 

{flfc, aj}2 = T^kjO-kaj (13.6) 

with the coefficients 

^fci = ^ {h,j+m - h+m,j + wSrnJ+m " w'f^j+m ' '^^klmj + ^kj ) (13.7) 

Here, in turn, the coefficients w^f!j are given in the A^-periodic case with g.c.d.(A^, p) = 1 by 
the formula 

k = j (mod A^) 

^S = \ " ' ' (13.8) 




k ~ j = nN (mod p), 1 < n < p — 1 



The Poisson structure (|13.6|) is nonlocal unless p = 1 (the case of BLl(m)), when it is 
given by the following brackets: 

O^k+l} = — afcOfc+i , ... {Ofc, Clfc+m} = —O^kO'k+m (13.9) 
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In particular, for m = 1 we obtain the quadratic Poisson bracket for the Volterra lattice. 
The corresponding Hamiltonians are: 

N 

H(a) = ^flfc for BLl(m) (13.10) 

k=l 

N 

H(a) = akak+i . . . ak+p-i for BL2(p) (13.11) 

k=l 

N 

H(a) = Y.^k'^kli---aklp-i for BL3(p) (13.12) 

k=l 

13.2 Lax representations 



The Lax representations of the Bogoyavlensky lattices fall into the class considered in ||K1 
and were also specified in |]B|. They have the form 

f=[T,B] (13.13) 

with the spectral-dependent (in the periodic case) matrices T = T(a, A), B = B{a,X). 
Actually, the natural phase space of all these Lax representations is the same algebra g as 
before. For BLl(m) the matrices T, B are given by 

N N 

T{a, X) = \Y. Ek+i,k + A--" ^ a^Ek^k+m (13.14) 

k=l k=l 

N N 

Bia, A) = 7r+(T'"+^(a, A)) = Y.i^k + a^-i + • • • + ak.m)Ek,k + A'"+^ ^ E,+rn+i,k (13.15) 

k=l k=l 

for BL2(p) they are given by 

N N 

Tia, A) = AP ^ Ek+p,k + A-^ (13.16) 

k=l k=l 

N 

B{a, A) = -7r_ (TP+^{a, A)) = -A"^"^ ^ afcOfc+i • • • afc+p^fe.fc+p+i (13.17) 

k=l 

and for BL3(j9) they are given by 

N N 

T{a, \) = yY. Ek+p,k + A ^ ak+iEk+i,k (13.18) 

k=l k=l 
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N 



B{a, A) = -n(T-^+\a, A) = -X'^^' ^ a,l,a,l, . . . a.^i^^.^^+P-i (13.19) 



k=l 



All three Lax representations above may be seen as having the form 

t= [t, ±7i±{d^{T))\ 
with appropriate Ad-invariant functions ip on namely, 

ip{T) = ^;^;^(tr(T™+^))^ for BLl(m) (13.20) 
cp{T) = ^(tr(T^+i)) for BL2(p) (13.21) 

cpiT) = L^(tr(T-^+i)) for BL3(p) (13.22) 

p — l\ /o 

(it is easy to see that the values of these functions in coordinates coincide with ( p.3.10| )- 
( |13.12| ), respectively). 

These Lax equations allow an r-matrix interpretation in which a certain quadratic 



bracket on g and its Dirac reduction to the manifold of Lax matrices are involved. It may be 
demonstrated that in general the Dirac correction to the Lax equations of motion vanishes 
for the Hamilton functions of the form (p{T) = (^tr(T")^^, where n = {pi + mi)ni with 
some rii E Z and m = mid, p = pid, d = g.c.d.{m,p). In particular, this explains the Lax 
equations for the Hamilton functions ( |13.20| )-( |13.22| ). 

13.3 Discretization of BLl 

We are now in a position to apply the recipe of Sect. |^ to the problem of discretizing the 
Bogoyavlensky lattices. This was done for the first time in |p8|| . Closely related results 
were obtained in |LHO|| , ||PN|| by different methods. We reproduce here the results of 



without proofs. As usual, the construction gives automatically for all discrete time systems 
(called hereafter dBLl, dBL2, dBL3, respectively) the invariant Poisson structure, the Lax 
representation, the integrals of motion, the interpolating Hamiltonian flows etc. The maps 
of all three families are nonlocal, but we demonstrate how to bring them to a local form by 
means of a suitable change of variables. Of course, the local forms of the maps dBLl-dBL3 
are Poisson with respect to the Poisson brackets on i3(a) which are the pull-backs of the 
corresponding brackets on B{a). However, we did not succeed in finding more or less nice 
formulas for such pull-backs (the only exceptions - the Volterra and the modified Volterra 
lattices). 

As dBLl(m) we take the map described by the discrete time Lax equation 

f = B-^TB+, = H+(j + /iT'"+^) (13.23) 
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Theorem 13.1 ||S8|| . The discrete time Lax equation ( p.3.23| ) is equivalent to the map a 
described by the equations 

5fc = ^^^ak (13.24) 
Pk 

where the functions j3k = /5fc(a) = 1 + 0{h) are uniquely defined for h small enough by the 
recurrent relations 

m / L \ 

h-hau = X{[l+ ^ (13.25) 
and have the following asymptotics: 

m 

Pk = l + hY,ak-j + Oih^) (13.26) 

j=0 



Remark. The factor is of the form 

N N 

B+ia, A) = n+(/ + /ir™+i) = J2 PkEk,k + /iA™+i J2 Ek+m+i,k (13.27) 

k=l k=l 

The map ( p.3.24| ) is nonlocal due to the functions f3k- In the simplest case m = 1 one has 
the continued fractions, terminating for the open-end boundary conditions: 

/3k-hak = l + — (13.28) 

■ + 



1 + hai 



and periodic for the periodic boundary conditions. For m > 1 there lacks even such an 
expressive mean as continued fractions to represent these nonlocal functions. 
The localizing change of variables B(a) ^ B{a) for dBLl(m) is: 

m 

ak = akWil + hAk-j) (13.29) 
i=i 

Obviously, the map ( |13.29| ) for h small enough is a local diffeomorphism. 

Theorem 13.2 The change of variables ( |13.29| ) conjugates the map dBLl(m) with the fol- 
lowing one: 

m m 

Sfc n (1 + hsik-j) = afc n (1 + hsik+j) (13.30) 
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Proof. Let us define the quantities (3k by the relation 



= l[{l + h8ik-j) 

j=0 



^13.31) 



and prove that they satisfy the recurrent relations ( |13.25[ ). Indeed, from ( |13.29| ), ( |13.31| ) it 
follows: 

Pk - hak = Y[{1 + hsik-j) = — 



Hence = ak/{l3k — hak), which, being substituted in the previous formula, implies ( |13.25| ) 



Now the uniqueness of solution of this latter recurrent system yields ( 13.31 ). Plugging 
( |I331| ), ( |13:29| ) into the equations of motion ( ^3:2^ ) results in ( p^ . ■ 



Theorem 13.3 The pull-back of equations of motion ( p.3.1| ) under the change of variables 



13.291 ) is given by the following formula: 

n,"li(l + hak+,) - n,"li(l + hsik-i, 



kk = SLk(l + hSLk 



h 



;i3.32) 



Proof. By a direct calculation one checks that the equations ( |13.32|) under the change of 
variables ( |13.29| ) are mapped to the equations of motion ( [13.11 ). ■ 



13.4 Discretization of BL2 

As dBL2(p) we take the map described by the discrete time Lax equation 

f = C.TCz\ C. =U_(l + hTP+^' 



;i3.33) 



Theorem 13.4 JS^. The discrete time Lax equation ( |13.33| ) is equivalent to the map a a 
described by the equations 

ak-h-fk-p /1QQ/l^ 
ak = ak+p+i (13.34) 

dk+p+l — illk+l 

where the functions 7^ = 7A:(a) = 0(1) are uniquely defined for h small enough by the 
recurrent relations 

ttk - h-fk-p = p (13.35) 



and have the asymptotics 



1 + hY[{ak-j - h-fk-p-j] 



j=0 



13.36) 
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Remark. The factor G_ is of the form 

N 

C.{a,X) = U^{l + hTP+')=I + hX-^P+'^Y.^kEk,k+p+i 

k=l 



;i3.37) 



The quantities 7^ render the equations of motion ( |13.34D nonlocal. The localizing change 
of variables for dBL2(p) reads: 



V i=i 

As usual, its bijectivity follows from the implicit functions theorem. 



[13.38) 



Theorem 13.5 The change of variables ( |13.38| ) conjugates the map dBL2(p) with the fol- 
lowing one: 

ak(^l + h f[ afc_, j =sik(^l + hf[ afc+, j (13.39) 

Proof. We proceed according to the by now already usual scheme. Define the quantities 7^ 
by the formula 



7fc = n ^fc+i 

j=0 



;i3.40) 



Then we immediately derive: 

flfe - h'jk-p = afc (13.41) 

and plugging this expression for a^ into ( |13.38D shows that the recurrent relations ( p.3.35| ) 
are satisfied. The uniqueness of the system of functions — hjk-p satisfying these relations 
justifies the expressions ( 13.40 ). Now putting ( 13.38[ ), ( |13.41 ) into the equations of motion 
( |13.34| ) allows to rewrite them as ( |13.39|) . ■ 



Theorem 13.6 The pull-back of the equations of motion ( |13.2|) under the change of vari- 
ables ( |13.38| ) is given by the following formula: 



afc = afc I ]J a-k+j - Yl ^k-j Yi ( ^ + ^ IT ^k-n+i 

\j=l j=l I n=l ^ i=l 



;i3.42) 



Proof. A direct, though tedious calculation shows that the equations of motion (|13.42|) are 
mapped on ( |13.2| ) by means of the map ( |13.38| ). ■ 
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13.5 Discretization of BL3 

As dBL3(p) we define the map with the discrete time Lax representation 

T = A_TAZ\ A_ = n_ (/ + hT-P+^) (13.43) 



Theorem 13.7 ||S8|| . The discrete time Lax equation ( [I3.43 ) is equivalent to the map a^a 

described by the equations 

ttfc - hak-p , 
au = r ak+p-i (13.44) 

where the functions ak = ak{a) = 0(1) are uniquely defined for h small enough by the 
recurrent relations ^ 

ak = n ]- (13.45) 

j^i (ik+j — nak+j^p 

and have the asymptotics 

p-i 

ak=ll<ji^ + 0{h)) (13.46) 
j=i 

Remark. The factor A_ is of the form 

N 

A_(a, A) = n_(/ + HT-P^') =1 + hX'P^' J2 ^kEk,k+p-i (13.47) 

fc=i 

The nonlocality of the equations of motion ( |13.44| ) is due to the functions ak{a). For 
example, for p = 2 they can be expressed as periodic continued fractions of the following 
structure: 

hak = 



h 

o-k+i 



h 

ak 

o-k-i — . _ 

The localizing change of variables for dBL3(p) is given by: 

afc = afc l^l + ZiHafc^ j (13.48) 
The bijectivity of this map is assured by the implicit functions theorem. 
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Theorem 13.8 The change of variables (|13.48|) conjugates the map dBL3(p) with the fol- 
lowing one: 



p-i 



p-i 



a J 1 + n a^^. = a J 1 + n a.^- 



;i3.49) 



j=0 



Proof. Defining the quantities ak by the formula 

p-i 

«fc = n ^kij 

we obtain with the help of ( |13.48D : 

Ofc - hak-p = a-k 



j=0 



;i3.50) 



;i3.51) 



Substituting this expression into (|13.50| ), we see that the recurrent relations ( |13.45|) are 
satisfied, which proves ( |13.50| ). Substituting ( |13.48| ), ( |13.51| ) into ( |13.44D , we immediately 
arrive at the equations of motion (p.3.49|) . ■ 



Theorem 13.9 The pull-back of the equations of motion ( |13.3|) under the change of vari- 
ables ( |13.48| ) is given by the following formula: 



fp-i 



p-i 



p-i 



n a^i - n ^k-j n 1 + ^ n ^k-n+i 



n=0 



p-l 



j=0 



^13.52) 



Proof - by direct (but tiresome) calculations. I 

We close the discussion of the local equations of motion for the dBL3(p) by noticing that 
under the change of variables a^ i-^ = a^^, h ^ —h the map ( |13.49| ) turns into 



qfc 1 1 - /i n ^fc-j 

3=0 



-1 



p-l 

qfc 1 1 - n ^fc+j 

3=0 



-1 



;i3.53) 



which is a local integrable discretization of the system ( |13.4| ), while the differential equations 
( |13.52| ) under the change of variables a^ i-^ = q^^, h ^ —h, t —t go into 



p-l 



n=0 



p-l 



(p-l p-l 
n ^k+j - n ^k-j I n ( 1 - ^ n 
j=i j=i 



■n+i 



1=0 



;i3.54) 



which is an integrable one-parameter deformation of ( |13.4| ). 
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13.6 Particular case p 



The case p = 2 of the Bogoyavlensky lattice BL3 is remarkable in several respects. The 
equations of motion for this case read: 



;i3.55) 



The localizing change of variables for the discretization of the system ( |13.55[ ) obtained in 
the previous subsection, is: 



afc = afc 1 + 



h 



afc + 



In the variables we have the following discretization: 

afc — a^ = 



;i3.56) 



;i3.57) 



(When considered on the lattice {t,k), the equation ( |13.57| ) is equivalent to the so called 
lattice KdV, which is a very popular object nowadays, cf. ||NC|| , ||NS|| , and references therein). 
The localizing change of variables ( 13.56 ) brings the system ( |13.55 ) itself into the form 



1 



1 



h 



,afc+i ^k~i, 
( 1 



afe-ia^^ 
1 



1 + 



h 



afc+iafc + /i Sik^k-i^h^ 



^13.58) 



The special properties of this case begin with the following observation. The subset of g 
consisting of the Lax matrices 

N N 

T{a, A) = ^ Ek+2,k + A XI ak+iEk+i,k (13.59) 

fc=l k=l 

is a Poisson submanifold for the corresponding quadratic r-matrix bracket, so that the Dirac 
reduction is not needed in giving an r-matrix interpretation to the bracket 

{ak,aj}2 = Trkjakttj (13.60) 

The coefficients TCkj are given (in the case of odd A^) by the formula ( |12.29| ). 
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Moreover, not only the quadratic r-matrix bracket, but also the linear one may be 
properly restricted to the set of the matrices ( p.3.59| ). The coordinate representation of the 
induced bracket on this set is given by: 



-1 



;i3.61) 



The two Poisson brackets ( |13.6U| ) and ( |13.61D are compatible, hence the system ( p.3.55| ) 
and its discretization given by Theorem |13.7| with p = 2 are bi-Hamiltonian. Obviously, the 
Hamilton function of the system (|13.55|) in the Poisson bracket (|13.61|) may be taken as 



Ho(a) 



N 
k=l 



;i3.62) 



It is easy to check that this function is a Casimir of the quadratic bracket. 

Being bi-Hamiltonian, the system (|13.55D and its discretization admit also an arbitrary 
linear combination of the brackets ( [L3.60| ) and ( |13.61|) as an invariant Poisson structure. A 
further remarkable feature is the following: there exists a linear combination of these two 
brackets whose pull-back under the map ( |13.56| ) allows a nice representation in the variables 



Theorem 13.10 The pull-back of the bracket 

{■,-}2 + M-,-}i 



on 13(a) under the map (|13.56|) is the following bracket on B{a) 



[13.63) 



[13.64) 



with the coefficients ( |12.29| ) . The map (|13.57| ) and the flow ( |13.58| ) are Poisson with respect 
to this bracket. 

Proof. A direct calculation shows that the brackets (|13.64|) for the variables a^ are sent by 
the map ( |13.56| ) into the brackets ( |13.63| ) for the variables a^. ■ 

This theorem allows also to derive the equations of motion ( |13.58| ) in a Hamiltonian 
manner. Indeed, these equations describe the Hamiltonian flow with the Hamilton function 

N 

h-^ log (afc(l + ha^^a^^)) 

k=l 

in the Poisson brackets (|13.64|) . Indeed, this function is a pull-back of Z^fcLi log(aA:), 
L3.55 p in the bracket (|13.63D. A direct calculation shows that this Hamil- 



which generates 

tonian system is governed by the differential equations (|13.58|) 

Remark. The results of this subsection agree with the results of the subsection |12.4 
after the change of variables qk = ^k^, i— 
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14 Alternative approach to Volterra lattice 



Starting from this point, we consider the systems with Lax representations in the direct 
products g = g^"^ rather than in g itself. We start with an alternative approach to the 
Volterra lattice based on a Lax representation in g ® g. 



14.1 Equations of motion and bi— Hamiltonian structure 

The version of VL we consider here is: 

Uk = Ukivk - ffc-i), Vk = Vkiuk+i - Uk) (14.1) 

Usually we let the subscript k change in the interval 1 < k < N and consider either open- 
end boundary conditions (fo = un+i = 0) or periodic ones (all indices are taken (mod A^). 
We consider mainly the periodic case, because the open-end one is similar and more simple. 
The relation to the form ( |10.1D is achieved by re-naming the variables according to 

Uk ^ a2k-i, Vk ^ 0'2k (14.2) 

So, in the present setting the A^-periodic VL consists of 2N particles. The phase space of 
VL in the case of periodic boundary conditions is the space 

W = M.^''{ui,vi,...,un,vn) (14.3) 

Two compatible local Poisson brackets on W invariant under the flow VL are given by 
the relations 

{ii'k,Vk}2 = -UkVk, {vk,Uk+i}2 = -VkUk+1 (14-4) 

and 



;i4.5) 



{Uk,Vk}3 = -UkVk{Uk + Vk), {Vk,Uk+l}3 = -VkUk+l{Vk + Uk+l) 

{Uk,Uk+l}3 = -UkVkUk+1, {Vk,Vk+l}3 = -VkUk+lVk+1 

respectively. The corresponding Hamilton functions for the flow VL are equal to 

N N 

}l,{u,v) = Y.Uk + Y.^k (14.6) 

k=l k=l 

and 

N N 

^o{u,v) = Y,^og{uk) or }lo{u,v) = "^log{vk) (14.7) 

fc=l k=l 

(the second function makes sense only in the periodic case; the difference of these two 
functions is a Casimir of {■, -js whose value is fixed to oo in the open-end case). 
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14.2 Lax representation 

Consider the following two matrices: 



N 



N 



N 



:i4.8) 



k=l 



k=l 



k=l 



These formulas define the "Lax matrix" {U, V) : W ^ g = g ^ g. 



Theorem 14.1 (see also The flow ( |14.1[ ) admits the following Lax representa- 

tion in g ® g: 



U = UC+ - B+U = B^U - UC^ 
V = VB+ - C+V = C-V - VB^ 



(14.9) 



with the matrices 



N N 

B+{u,v,\) = '^{uk + Vk-i)Ekk + X'^Ek+i^k 

k=l k=l 
N N 

C+{u,v,\) = ^{uk + Vk)Ekk + \J2^k+i,k 

k=l k=l 

N 

5_ (u, v,\) = X'^Yj UkVkEk,k+i 

k=l 
N 

C- (m, f , A) = A~^ X Uk+iVkEk,k+i 

k=l 

Corollary. The matrices 

T+{u, V, X) = U {u, V, X)V{u, V, A), T^{u, v, A) = V{u, v, X)U{u, v, A) 
satisfy the usual Lax equations in g: 



(14.10) 
(14.11) 
(14.12) 
(14.13) 

(14.14) 
(14.15) 



The matrices T± are easy to calculate explicitly. From the corresponding formulas one 
sees that the matrices B±, C± allow the following representations: 



B, 



TT 



+ (T+) = 7r+(f/y), B_ = 7r+(T_ 



7r_ 



-{VU) 



C+ = 7r_(T+) = 7i^{UV), = 7r_(r_) = n^{VU) 

The Lax equations ( |14.9| ) may be given an r-matrix interpretation ||S12|| in terms of a 
certain quadratic bracket on g = gr ® gr, the induced bracket on W being {-, ■}2. 
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14.3 Discretization 

To find an integrable time discretization for tlie flow VL, we apply the recipe of Sect. |]witli 
F{T) = I + hT, i.e. we consider the map described by the discrete time "Lax triads" 

U = B-^UC+ = B^UCZ} , V = C-^UB+ = C_VB_^ (14.16) 

with 

B± = n± (/ + huv) , c± = n± (/ + hvu) 



Theorem 14.2 The discrete time Lax equations ( |14.16| ) are equivalent to the map {u,v) t— > 
(m, v) described by the following equations: 

Uk = Uk—, Vk = Vk (14.17) 

Pk Ik 

where the functions (3k = Pk{u, v) = 1 + 0{h) , •jk = 1k{u, v) = 1 + 0{h) are uniquely defined 
for h small enough by the recurrent relations 

Pk - huk = ^'-l = 1 + (14.18) 

Ik - hvk = — = 1 + — 14.19 

(Jk - huk Pk - huk 



and have the asymptotics 



(5k = l + h{uk + Vk-i)+0{h^) (14.20) 
7fc = l + h{uk + Vk)+0{h^) (14.21) 



Remark. The matrices B±, C± have the following expressions: 

AT TV 



B+{u,v,\) = U+[l + hUV) = PkEk,k + hXY, Ek+i,k (14.22) 

k=l k=l 
N N 

C+{u,v,\) = U+[l + hVU) = Y.lkEk,k + h\Y.Ek+i,k (14.23) 



k=l k=l 

N 



B4u,v,\) = U_[l + hUV) = I + h\-'Y.^Ek,k+i (14.24) 

k=i Pk 



N 

C4u,v,X) = U^[l + hVU) = I + h\-^Y.^^^^Ek,k+i (14.25) 

k=i 
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Proof. The general structure of the factors B±, C±, as given in ( |14.22|) -( p.4.25| ), as well 
as the expressions for the entries of -B_, C_, follow directly from the defining equalities 
-B+S_ = / + hUV, C+C- = / + hVU. For the entries (5k, ■jk of C+ one obtains the 
following recurrent relations: 



pk = l + h{uk + Vk-i) 



h'^Uk-iVk- 



7fc = 1 + h{uk + Vk) - 



h'^UkVk-i 



Pk-i 1k-i 
Now notice that these relations coincide with (|10.9|) after re-naming ( |14.2|) and 



(14.26) 



Hence we may use the proof of Theorem |10.1| to establish the alternative recurrent relations 
( |141^ ), ( ImTqI) . The equations of motion ([liT^) follow directly from B+U = UC+, C+V = 
VB^. It is important to notice that ( |14.17 ) become identical with ( 10.11|) after the above 
mentioned re-namings of variables. This allows to denote consistently the map constructed 
in this theorem by dVL. ■ 



14.4 Local equations of motion for dVL 

Now we can simply reformulate the results of Sect. |TU|in our new notations. 

The localizing change of variables W(u, v) i— VV('U, v) for dVL is given by the following 
formulas: 

Mfc = Ufe(l + /ivfc_i), Vk = ^k{l + hu.k) (14.27) 

Due to the implicit function theorem, these formulas define a local diffeomorphism for h 
small enough. 

Theorem 14.3 The change of variables ( |14.27| ) conjugates the map dVL with the map 
(u, v) ^ (u, v) governed by the following local equations of motion: 

Ufc(l + /iVfc_i) = Ufc(l + /iVfc), Vfc(l + /lUfc) = Vfc(l + /lUfe+i) (14.28) 



Let us mention that the functions l3k{u,v), 'jk{u,v) in the localizing variables are given 
by the formulas 

A = (l + /iufe)(l + /ivfc_i), 7fc= (l + /ivfc)(l + /iufc) (14.29) 

so that 

(3k - huk = 1 + hvk-i, -Ik - hvk = I + huk (14.30) 
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Remark. Cosidering the equations ( 14.28| ) as lattice equations on the lattice {t,k), and 
performing a linear change of independent variables, one arrives at the explicit version of 
dVL, cf. [Vll, iFg. 



Theorem 14.4 The pull-back of the bracket 

{•,-}2 + M-,-}3 (14.31) 
on }V{u,v) under the change of variables ( |14.27| ) is the following bracket on >V(u, v); 
{ufc,Vfc} = -UfcVfc(l+/iUfc)(l+/iVfc), {vfc,Ufc+i} = -VfeUfe+i(l+/iVfc)(l+/iUfe+i) (14.32) 
The map (|14.28|) is Poisson with respect to the bracket (|14.32| ). 



Theorem 14.5 The pull-back of the flow VL under the map ( |14.27| ) is described by the 
following equations of motion: 



ilk = Ufe(l + /iUfe)(vfe - Vfc_ 



Vfe = Vfc(l + /lVfe)(Ufc+i - Vfc) 



14.33) 



14.5 Lax representation for VL2 

We consider now the g^g formulation of the second flow VL2 of the Volterra hierarchy. We 
use the notations Pk, Ik, B±, C± etc. for objects analogous to those relevant for VL without 
danger of confusion. 

The flow VL2 is described by the following differential equations on W: 

Uk = Uk{vk{Uk+l + Vk + Uk) - Vk-l{Uk + Vk-l + Uk-l)) 

Vk = Vk{uk+i{vk+i + Uk+i + Vk) - Uk{vk + Uk + Vk,)) (14.34) 
The Hamilton functions of this flow are: 

N N 

H2 (m, W) = - X! + ^fc) + H {Uk+lVk + VkUk) 
^ k=\ k=l 

in the quadratic bracket {■, ■}2; and B.i{u,v) in the cubic bracket {■, - ja. 

The Lax representation for the flow VL2 is of the type (|3.2| ) with m = 2 and /(T) = T^. 



54 



Theorem 14.6 ||K1|| . The flow ( p.4.34| ) admits the following Lax representation in g ® g: 

U = UC+ - B+U = B^U - UC^ 

(14.35) 

V = VB, - C,V = C^V - VB^ 



with the matrices 

N 

B+{u,v,X) = J2{(^k + Vk-i)'^ + UkVk + Uk-iVk-i)Ekk 

k=l 

N N 

+ A ^(Ufc+l + Vk+Uk + Vk-l)Ek+l,k + X! Ek+2,k 
k=l k=l 

N 

C+{u,v,X) = ^{^{uk + Vk)'^ + Uk+iVk + UkVk-i)Ekk 

k=l 

N N 

+A ^(wfc+i + Vk+i +Uk + Vk)Ek+i,k + ^ Ek+2,k 

k=l k=l 

N 

B_ (u, V, A) = A~^ J2 ("fc+i + '"k + Uk + Vk-i)ukVkEk,k+i 

k=l 
N 

+A~^ ^ Uk+iVk+iUkVkEk,k+2 

k=l 

N 

C-(u, V, A) = A~^ "^{uk+i + Vk+i + Uk + Vk)uk+iVkEk,k+i 

k=l 

N 

+ A'^ ^ Uk+2Vk+lUk+lVkEk,k+2 

k=l 



Obviously, the expressions for B^, (7+ may be obtained from ( |9.4| ) with the help of the 
substitutions (|14.47|) , (|14.48|) , respectively, and C_ follow analogously from (|9.5| ). 



14.6 Discretization of VL2 

Applying the recipe of Sect. ^ with F{T) = I + hT"^, we take as a discretization of the flow 
VL2 the map described by the discrete time Lax triads 

U = B-^UC+ = B^UCZ^ , V = CfUB+ = C VBZ^ (14.36) 

with 

B± = U^[l + h{UVy), C± = Il±[l + h{VUf 
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Theorem 14.7 The discrete time Lax equations ( |14.36| ) are equivalent to the map {u,v) i— > 
(tt, v) described by the following equations: 

Uk = Uk — , Vk = Vk (14.37) 

Pk Ik 

where the functions Pk = Pk{u, v) = 1 + 0{h) , •jk = 1k{u, v) = 1 + 0{h) are uniquely defined 
for h small enough simultaneously with the functions Sk = 0(1), = 0(1) by the system of 
recurrent relations 

Pk - h{6k - Uk+i)uk = \ (14.38) 

7yt_i - h{ek-i - Vk)Vk-i 

Ik - h{ek - Vk+i)vk = "5 TT7^ ^ — (14.39) 

Pk - h{dk -Uk+i)uk 

A II, huk-iVk-i5k-i (^AA(^\ 
bk = Uk+i +Vk + Uk + Vk-i (14.40) 

Pk-l 

hukVk-itk-i (^AA^\ 
ek = Uk+i + Vk+i +Uk + Vk (14.41) 

7fe-i 

The auxiliary functions (3k, jk have the asymptotics 

Pk = 1 + h(^{uk + Vk-iY + UkVk + Uk^iVk-i) + 0{h^) (14.42) 
7fc = 1 + h(^{uk + VkY + Uk+iVk + UkVk-ij + 0{h'^) (14.43) 
Remark. The matrices B^, have the foUowing expressions: 

N N N 

^k+2,k 



B+{u,v,\) = U+(l + h{UVY) = J2PkEk,k + h\Y,SkEk+i,k + h\^Y.E. 

k=l k=l k=l 

N N N 

C+iu,v,\) = n+(l + hiVUf) = Y.lkEk,k + hXY,ekEk+i,k + hX^Y.E. 



k+2,k 

k=l k=l k=l 



Proof. The scheme of the proof is standard. First of aU, the general structure of the factors 
B±, C± is clear from the bi-diagonal structure of the matrices U, V. The recurrent relations 
for the entries of the matrices B+, C+ follow, in principle, from ( |9.8|) , ( |9.9|) with the help 
of substitutions ( 14.47| ), ( |14.48 ). It is easy to see that these relations coincide with ( 11.8 



( |11.6| ) after re-naming (|14.2|) and 

Pk ^ /?2fc-l, Ik ^ /?2fc, h ^ 62k-l, Cfe S2k 



The equations of motion ( |14.37|) , following from B^U = C+V = VB+, also coincide 



with ( 11. 4| ) after the above re-namings. Thus the discretization introduced in this theorem 



agrees with the one from Sect. 11 and may be consistently denoted by dVL2. 
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14.7 Local equations of motion for dVL2 



Here we translate the corresponding results from Sect. |TI| into our present notations. The 



localizing change of variables for the map dVL2 is given by the formulas 
(l + /ivti) _ (1 + hul) 

(1 - /lUfcVA._i)(l - hVk-lUk-l) ' (1 - /lVfeUA.)(l - hUkVk-l 



Theorem 14.8 The change of variables ( |14.44|) conjugates the map dVL2 with the map 
(u, v) (u, v) described by the following local equations of motion: 

(1 + /^vLi) ^ {l + h^rD 

Ufc J- , ~ ~ rrz 7^ ~ 7 — 



1 - /iUfcVfc_i)(l - hvk-iUk-i) (1 - /iUfc+iVfc)(l - hvkUk) 

(14.45) 



;i + hul) (1 + huj 



Vfc , 7-^-3 T = Vfc 



1 - hvkUk){l - hukVk-i) (1 - /iVfe+iUfc+i)(l - /lUfc+iVfc) 



Theorem 14.9 The pull-back of the flow VL under the change of variables ( |14.44|) is de- 
scribed by the following equations of motion: 



lifc = Uk{l + hul) 



Vfc = Vfc(l+/lV^) 



1 - hvkUk 1 - hukVk-i J 

(14.46) 



1 — huk+i^'k 1 — hvkUk 



14.8 Miura relations to the Toda hierarchy 

We have seen in Sect. ^ that the flow VL may be considered as a restriction of the second 
flow TL2 of the Toda hierarchy. There exists a relation of a completely different nature with 
the flow TL. Namely, the flow VL is Miura related to the first flow TL of the Toda hierarchy, 
while the flow VL2 is Miura related to TL2. To see this notice that the matrices T±{u, v, A) 
from (|14.14D have the same tri-diagonal structure as the Toda lattice Lax matrix ( ^.9| ) with 



ak = UkVk, bk=Uk + Vk-i (14.47) 



or 



ak = Uk+iVk, bk=Uk + Vk (14.48) 

respectively. These two pair of formulas may be considered as two Miura maps A4± : 
W{u,v) I— i> T{a,b). The following holds (see 
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1. The both maps A4± are Poisson, if W is equipped with the bracket ( |14.4|) , and T is 
equipped with the bracket 



2. The both maps A4± are Poisson, if W is equipped with the bracket ( |14.5|) , and T is 
equipped with the bracket 



3. The flow VL (|T4T|) is conjugated with the flow TL and the flow VL2 ([HIS^ ) is 
conjugated with the flow TL2 (|9.1|) by either of the maps A4±. 

We now translate these statements into the language of localizing variables. Since we have 
two different localizing changes of variables (for dVL and dVL2), two different translations 
are necessary. 

We start with the case of the localizing change of variables ( [14.27| ). 



Theorem 14.10 [^| 

a) Define two maps M.± : W(u, v) 

and 



T(a, b) by the formulas 

l + hhk = {1 + huk){l + hvk-i) 



afc = Ufc+iVfc, 1 + hhk = {1 + huk){l + hvk) 
respectively. Then the following diagram is commutative: 



(14.49) 
(14.50) 



>V(u,v) 



Tfa,b) 



( imri) 



(Km 



W{u,v) 



M 



± 



T{a,b) 



b) The both maps M.± are Poisson, z/ >V(u, v) is equipped with the bracket ( |14.32| ), and 
T(a, b) is equipped with the bracket (|8.25| ). 
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c) The local form of the dVL ( |14.28| ) is conjugated with the local form of the dTL (|8.2CI| ) 
by either of the maps M-|-. 



Proof. The first statement is verified by a direct check, the second and the third ones are 
its consequences. ■ 



In |[K1|| this theorem was formulated without any relation to the problem of integrable 
discretization. In the context of discrete time systems the formulas ( |14.49| ), ( |14.5U| ) were 
found also in ||HT|] , however, without discussing Poisson properties of these maps. 

Concerning the localizing change of variables ( |14.44| ) for dVL2, we come to the following 
results. 

Theorem 14.11 Define the maps Nj. : VV(u, v) T(a, b) by the formulas 



and 



Ufc+lVfc Ufc + Vfc 

Sik — ; ; Ofc — 

1 - /lUfc+iVfc 1 - /lUfcVfc 

respectively. Then the following diagram is commutative: 



:i4.5i) 



;i4.52) 



>V(u,v) 



T(a, b) 



T44|) 



T{a,b) 



M- 



and the local form of the dVL2 (|14.45|) is conjugated with the local form of the dTL2 (|9.15| ) 
by either of the maps Nj. . 

Proof - by a direct check. ■ 
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15 Relativistic Toda lattice 

15.1 Equations of motion and tri— Hamiltonian structure 

The relativistic Toda lattice was invented by Ruijsenaars ||Ruij|| , and further studied in 



BK1|| , ||(J.b'ZK|| and numerous other papers. In particular, the tri-Hamiltonian structure was 



elaborated in the latter reference. 

We consider here two flows of the relativistic Toda hierarchy: 

dk = dk{ck - Ck-i), Ck = Ck{dk+i + Ck+i - dk - Ck-i) (15.1) 

and 

4 = 41 -r^ 1 , Ck = Ck{-r- -1 1 (15.2) 

The both systems may be considered either under open-end boundary conditions {cq = cn = 
0), or under periodic ones (all the subscripts are taken (mod N), so that d^+i = di, 4 = djy, 
Civ+i = Ci, Co = Cn). We shall denote the first flow by RTL+, and the second one by RTL— . 
The phase space of the flows RTL± in the case of the periodic boundary conditions may be 
defined as 

n = R^^{di,Ci,. . . ,dN,CN) (15.3) 

This space carries three compatible local Poisson bracket, with respect to which the flows 
RTL± are Hamiltonian. 

The linear Poisson structure on TZ is defined as 

{4,Cfc}i = -Cfc, {cfc,4+i}i = -Cfc, {4,4+i}i = Cfc (15.4) 

The quadratic invariant Poisson structure on TZ is given by the brackets 

{dk,Ck}2 = —dkCk {Cfc, 4+1)2 = —Cfc4+1, {Ck,Ck+l}2 = —CkCk+l (15.5) 

Finally, the cubic Poisson bracket on TZ is given by the relations 

{4,Cfc}3 = — 4cfc(4 + Cfe), {cfc, 4+1)3 = — Cfc4+i(cfc + 4+1), 

{4,4+1)3 = —dkCkdk+i, {cfe,Cfc+i)3 = — CfcCfe+i(cfc + 24+1 + Cfc+i), 

{dk,Ck+i}3 = —dkCkCk+i, {cfc, 4+2)3 = —CkCk+idk+2, 



{Cfc, 0^+2)3 — —CkCk+lCk+2 



;i5.6) 
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The Hamilton functions for the flow RTL+ in the brackets (|15.4| ), (|15.5|) , ( |15.6|) are equal to 
H2 (c,(i), Hi (c,(i), and B.o{c,d), respectively, where 

N N 

Mc,d) = -E(4 + Cfc)' + E(4 + Cfc)cfc_i (15.7) 

^ fc=l k=l 

N 

Hi(c,d) = E(4 + Cfc) (15.8) 

k=l 
N 

Ho(c,ci) = Elog(^fc) (15.9) 

k=l 

Similarly, the Hamilton functions for the flow RTL— in these three brackets are equal to 
— Ho(c, (i), H_i(c, (i), H_2(c, (i), respectively, where 

TV r _|_ 

H-i(c,rf) = (15.10) 

w ( ^\ 1 ^ (4 + Cfc)^ , ^ (4-1 + Ck-i)ck 

H_2(c,rf) = ^E ^2 ,2 ^ 323 (15.11) 

^ A;=l "A:+1"A: fc=l "A:+l«fc«fc-l 

15.2 Lax representation 

The most natural Lax representation for the relativistic Toda hierarchy is the one living 
in g = g ^ g ||S3|| , which is in many respects analogous to the Lax representation for the 
Volterra hierarchy considered in the previous section. 
Introduce the matrices 

TV TV TV 

U{c, rf. A) = E dkEkk + A E Ek+i,k, V{c, d,X)=I- A-i E CkEk,k+i (15.12) 

k=l k=l k=l 

Theorem 15.1 | ]S3| j . The equations of motion ( |15.1D are equivalent to the following Lax 
equations in g ® g (Lax triads): 

tl = UB-AU, V = VB-AV (15.13) 

where 

TV TV 

A(c,rf,A) = E(^fe + Cfc-i)^fcfc + AE^fc+i,fc (15.14) 

fc=l k=l 

TV TV 

B{c,d,\) = Y.(dk + Ck)Ekk + \Y.Ek+i,k (15.15) 

k=l k=l 
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The equations of motion (|15.2|) are equivalent to the following Lax triads: 

il = UD-CU, V = VD-CV 

where 



N 



C{c, d, A) 
D{c, d,\) 



Ck 



k=i ^k+i 



■E, 



N 



k=l 



dk 



k,k+l 



k+l 



Corollary. The matrices 

T+{c,d,X) = U{c,d,X)V-\c,d,X), T_{c,d,X) = V-\c,d, X)U{c,d, X) 
satisfy the usual Lax equations in g. Namely, for the flow ( |15.1| ) 

f+ = [T+,A], f_ = [T_,B] 

and for the flow ( |15.2D 

t+ = [T+,C], t_ = [T_,D] 
The following formulas are easy to establish by a direct calculation: 

A = n+iT+), B = n+{T^) 
C = 7r_(T-i), D = n4T:') 
Hence the Lax equations (|15.20|) and ( |15.21| ) may be presented as 



(15.16) 

(15.17) 
(15.18) 

(15.19) 

(15.20) 
(15.21) 

(15.22) 
(15.23) 



T 



T,MT) 



and 



T 



T,7r_(T-i) 



respectively. These Lax equations in g may be given an r-matrix interpretation in the cases 
of the linear and of the quadratic Poisson brackets; for the Lax triads in gr (g) ^ an r-matrix 
interpretation is known in the case of the quadratic bracket only 



15.3 Discretization of the relativistic Toda hierarchy 

We see that actually the pairs (f/, V~^) satisfy the Lax equations of the form (|3.2|) (with 
m = 2). This allows to apply our general recipe to find integrable discretizations of the flows 
RTL±. This results in considering the following discrete time Lax triads: 

u = n;i(F(T+))-f/-n+(F(T_)) = n_(F(T+))-f/-n:i(F(T_)) 

V = n;i(F(r+))-\/.n+(F(r_)) = n_(F(T+))-v.n:i(F(T_)) 

with 

F(T) =I + hT and F(T) = / - HT'^ 

respectively. It turns out that for the flow RTL+ the version with the 11+ factors is more 
suitable, while for the RTL— flow the version with the n_ factors is preferable. 
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15.4 Discretization of the flow RTL+ 

Consider the discrete time Lax triad 

U = A-^UB, V = A-^VB (15.24) 

with 

A = n+(/ + /iT+), B = U+{I + hT^) 
implying also either of the two equivalent forms of a convenient Lax equation: 

f+ = A-^T+A, f. = B-^T^B (15.25) 



Theorem 15.2 ||S7|| . The equations (|15.24|) are equivalent to the map {c,d) i— > {c,d) de- 
scribed by the following equations: 

dk = dk — , Cfc = Cfc (15.26) 

where the functions = Ofc(c, d) = 1 + 0{h) are uniquely defined by the recurrent relation 

ak = l + hdk + (15.27) 

and the coefficients bk = bfc(c, d) = 1 + 0{h) are given by 

ak + hck ak+i — hdk+i , . 

bk = ak-i — = afc — (15.28) 

Ofc-i + iiCk-i Ofc — tidk 

The following asymptotics hold: 

ak = l + h{dk + Ck-i) + 0{h^) (15.29) 

bk = l + h{dk + Ck) + 0{h^) (15.30) 

Remark. The auxiliary matrices A, B are bi-diagonal: 

N N 

A(c, d, A) = n+(/ + hT+) = ^ akEkk + hXj2 Ek+i,k (15.31) 

k=l k=l 
N N 

B{c, d, A) = n+(J + hT_) = ^ bkE,, + h\Y^ Ek+i,k (15.32) 

k=l k=l 
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Proof. The general bi-diagonal structure of the factors A, B follows from their definition. 
The simplest way to derive the recurrent relations ( |15.27| ) for the entries of A is to notice 
that 

A = U+(l + hUV-^) = n+(V + hU), because V E 



The equations of motion ( |15.26| ) and the relations ( p. 5 .281 ) follow now from AU = UB, 

Av = VB. m 



The map ( |15.26| ) will be denoted dRTL+. Due to the asymptotic relations (|15.29|) 



( p.5.30| ) it clearly approximates the fiow RTL+. As usual, it is tri-Poisson, allows the same 
integrals and the same Lax matrix as the fiow RTL, but has a drawback of nonlocality. Its 
source are the functions a^, which in the open-end case may be expressed as terminating 
continued fractions: 

hck-i 

ak = l + hdk + 



1 + hdk-i + 



hci 



1 + hdi 

In the periodic case the a^'s may be expressed as infinite periodic continued fractions of an 
analogous structure. 

15.5 Local equations of motion for dRTL+ 

Introduce another copy of the phase space TZ parametrized by the variables c^, and 
consider the change of variables TZ{c, d) 7?.(c, d) defined by the following formulas: 

4 = dfe(l + /iCfc_i), Ck = Ck{l + hdk) {1 + hck-i) (15.33) 

Obviously, due to the implicit function theorem, for h small enough this map is locally a 
diffeomorphism. 



Theorem 15.3 The change of variables (|15.33|) conjugates dRTL+ with the map described 
by the following local equations of motion: 

dfc(l + hCk-i) = dfc(l + hCk) 

(15.34) 

Cfc(l + hdk){l + hCk-i) = Cfe(l + hdk+i){l + hck+i) 
Proof. The crucial point in the proof of this theorem is the following observation: the 



parametrization of the variables (c, d) according to (|15.33| ) allows to find the coefiicients 



(defined by the recurrent relations ( |15.27| )) in the closed form, namely: 

ak = il + hdk)il + hCk-i) (15.35) 
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Indeed, if we accept the last formula as the definition of the quantities a^, then we obtain 
successively from ( p.5.35[ ) and ( |15.33| ): 



afc = 1 + hdk{l + hck-i) + /iC/j_i = 1 + hd^ + 



hck- 



So, the quantities defined by ( |15.35| ) satisfy the recurrent relation ( |15.27| ), and due to the 
uniqueness of solution our assertion is demonstrated. Now (|15.28|) and ( |15.35|) yield: 



[1 + hdk){l + hck) 



(15.36) 



and the equations of motion (|15.34|) follow directly from ( p.5.26| ), ( |15.35| ), (|15.36| ). ■ 

It turns out that the pull-back of either of the brackets ( |15.4|) , ( |15.5|) , ( |15.6|) is highly 
nonlocal. Nevertheless, there exist certain linear combinations thereof, whose pull-backs are 
local. 



Theorem 15.4 a) The pull-back of the bracket 

{■,-}i + Mt}2 



;i5.37) 



on 7l{c, d) under the change of variables ( |15.33| ) is the following Poisson bracket on 
7^(c,d).■ 



{dfc,Cfc} = -Cfe(l + Mfe) 
{cfc,dfe+i} = -CA.(l + Mfe+i) 

I + hdk){l + hdk+i) 



^15.38) 



{dfc^dfc+i} — Cfc 



1 + hCk) 



b) The pull-back of the bracket 

{•,-}2 + Mt}3 (15.39) 

on TZ{c, d) under the change of variables ( p.5.33| ) is the following Poisson bracket on 
7^(c,d).■ 



{dfc,Cfc} = -dfcCfc(l + /idfc)(l + /iCfc) 
{Cfc,dfc+i} = -Cfcdfc+i(l + /iCfc)(l + /idfc+i) 
{cfc,Cfc+i} = -CkCk+i{l + hck){l + hdk+i){l + hck+i] 



(15.40) 



c) The brackets ( |15.38|) , ( |15.40|) are compatible. The map ( |15.34| ) is Poisson with respect 
to both of them. 
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Proof - by a direct check. Notice a curious feature of the bracket (|15.38|) : it is nonpolynomial 
in coordinates (though still local). ■ 



Theorem 15.5 The pull-back of the flow RTL+ under the map ( |15.33| ) is described by the 
following equations of motion: 

dfc = dfe(l + hdk){ck - Ck-i) 

(15.41) 

Cfc = Cfc(l + /iCfc)(dfc+i + Cfc+i + /idfc+iCfc+i - dfc - Cfc_i - /idfcCfc_i) 

Proof. We will use by the proof only the bracket of the part b) of the previous theorem. 
Obviously, the pull-back we are looking for is a Hamiltonian system with the Hamilton 
function, which is a pull-back of /i~^Ho(c, d) (indeed, this function is a Casimir function for 
{■, ■}2 and is the Hamilton function of the flow RTL+ in the bracket -ja). Calculating 
the equations of motion generated by the Hamilton function 

N 



h ^^log(dfc(l + hCk_ 



k=l 



in the Poisson brackets ( |15.40| ), we obtain ( |15.41| ). 



15.6 Discretization of the flow RTL— 

Consider the discrete time Lax triad 

U = CUD-\ V = CVD-^ (15.42) 

with 

C = H_(/-/iT-^), D = U4l-hT-^) 
implying also the convenient Lax equations: 

f+ = CT+C~\ f^ = DT^D-^ (15.43) 



Theorem 15.6 ||^. The equations (|15.42|) are equivalent to the map {c,d) i— > {c,d) de- 
scribed by the following equations: 



dk = d 



k+l 



Cfc 



(15.44) 
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where the functions t)^ = Dk{c,d) = 0(1) are uniquely defined by the recurrent relation 

= r^^T (15.45) 

dk — 11 — nUk-i 

and the coefficients = Ck{c,d) = 0(1) are given by 

dk - hdk-i Ck + hi3k i^r. , „x 

= '^^ 1 r~ = ^^=+1 TT (15.46) 

The following asymptotics hold: 

^k = ^ + 0{h), ck = -p^ + 0{h) (15.47) 

Remark. The auxiliary matrices C, D in the discrete time Lax equations admit the 
following expressions: 

AT 

Cic,d,X) = I + h\-'Y.'kEk,k+i (15.48) 

k=l 

N 

D{c, d,X)=I + hX-' ^kEk,k+i (15.49) 

fc=i 



Proof - analogous to that of Theorem 15.2 



The map ( |15.44| ) will be called hereafter dRTL— . Like dRTL+, it is tri-Poisson etc. It 
is nonlocal because of the presence of the functions 0^, which in the open-end case have the 
following finite continued fractions expressions: 



Cfc 

ffc = 



, , hck-i 
dk - h 



4-1 - h 



hci 



di — h 



15.7 Local equations of motion for dRTL— 

To bring the map dRTL— to the local form, another change of variables 7^(c, d) i-^ 7^(c, d) 
is necessary: 

Again, for h small enough this map is locally a diffeomorphism, due to the implicit function 
theorem. 
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Theorem 15.7 The change of variables (|15.50|) conjugates dRTL— with the map described 
by the following local equations of motion: 



1 + 



hci 



fl5.511 



Cfc 1 



Proof. This time the crucial component of the proof is the following remarkably simple local 
formula for the coefficients 0^ (defined by the recurrent relations ( |15.45| )) in the coordinates 
(cfe, dfe): 

ffc = ^ (15.52) 
dfc 

Indeed, if we use both ( p.5.50| ) and ( |15.52[ ) as definitions, then we obtain: 

h = dk — h — h^k-i 



Ck J 

— = dfc 



Hence, the quantities defined by ( |15.52| ) satisfy the recurrent relation ( |15.45| ), and due to 
the uniqueness of solution our assertion is proved. From ( p.5.46| ) and ( |15.52| ) we obtain also: 



Ck 



Ck 



dfc 



+1 



;i5.53) 



Now the equations of motion ( |15.51| ) follow directly from ( |15.44| ), ( |15.52| ), ( |15.53| ). ■ 

The Poisson properties of the change of variables (|15.50 ) are similar to that of ( 15.33 ). 
Namely, the pull-backs of either of the brackets ( |15.4|) , ( |15.5| ), (|15.6| ) are nonlocal, but there 
exist linear combinations thereof, whose pull-backs are local. 

Theorem 15.8 a) The pull-back of the Poisson bracket 

{■,-}2-M-,-}i (15.54) 
on TZ{c,d) under the change of variables ( |15.5CI| ) is the following bracket on Tl{c,d): 



{dfc,Cfc} 



-dfcCfc 1 - 



h 



{cfc, dfc+i} = -Cfcdfc+i 1 - 



h 



{Cfc,Cfc+i} = -CfcCfc+1 1 



h 



lfe+1. 



(15.55) 
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b) The pull-back of the Poisson bracket 

{■,-}3-M-,-}2 (15.56) 
on TZ{c,d) under the change of variables ( |15.5U| ) is the following bracket on TZ{c,d): 



{dfe,Cfe} = -dfcCfc(dfc + Cfc) ^1 - -^^ 
{cfc, dfc+i} = -Cfcdfc+i(cfc + dfc+i) ^1 - -J^ 
{dfc,dfc+i} = -dfcCfcdfc+i ( 1 - -J- ) ( 1 



{cfc, Cfc+i} = -CfcCfc+i(cfc + 2dfc+i + Cfc+i) ( 1 - I (15.57) 



dfc+i 



{dfc,Cfc+i} = -dfcCfcCfc+i ( 1 - ] ( 1 



dfc / \ dfc+i 



{cfc, dfc+2} = -CfeCfc+idfe+2 1 - :3 — 1 - :^ — 

{cfc, Cfc+2} = -CfcCfc+iCfc+2 (1 - 1 \ T~~ 



c) r/ie brackets (|15.55|) , (|15.57|) are compatible. The map ( |15.51| ) zs Poisson with respect 



to both of them. 
Proof consists of straightforward calculations. 



Theorem 15.9 The pull-back of the flow RTL— under the map (|15.5CI|) is described by the 
following equations of motion: 



dfcdfc+i \ dfcdfc_|_i/ dfc_idfc V dfc_id 



(15.58) 



Cfc = Cfc ; 1 



dfc dfc+i / V dfcdfe+i 



Proof. We use the part a) of the previous Theorem. The flow under consideration is a 
Hamiltonian system with the Hamilton function, which is a pull-back of h~^B.o{c, d) (indeed, 
this function is a Casimir function for {■, ■}2 and is a Hamilton function of the flow RTL— 
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in the bracket —h{-,-}i). Calculating the equations of motion generated by the Hamilton 
function 

hck \ 



N N / 

/i-i^log(dfc) + /i-i^log 1 + 

fc=l k=l \ 

in the Poisson brackets ( |15.55| ), we arrive at (|15.58| ). ■ 



15.8 Third appearance of the Volterra lattice 

It is interesting to remark that the flow RTL+ allows the reduction dk = 0, in which it turns 
into the Volterra lattice 

The Lax representation of the RTL+ flow survives this reduction, delivering a new (third) 
Lax representation for the VL. Also the quadratic and the cubic Poisson brackets ( |15.5| ) and 
( |15.6| ) allow this reduction and turn into the corresponding objects for the VL. 

It may be verified that the map dRTL+ in the reduction = turns into dVL, although 
this discretization is based on a completely different Lax representation. Naturally, the same 
holds for the local forms of these maps. 



16 Belov— Chaltikian lattice 

16.1 Equations of motion and bi— Hamiltonian structure 



In their studies of the lattice analogs of l^-algebras, Belov and Chaltikian PQ| , ||ABQ] 
found an interesting integrable lattice (hereafter BCL): 

bk = bk{bk+i - bk-i) - Ck + Ck-i, Ck = Ckibk+2 - bk-i) (16.1) 

This system may be viewed as an extension of the Volterra lattice (which appears as a = 
reduction of the above system). Belov and Chaltikian established also the bi-Hamiltonian 
structure of this system. Namely, its phase space, which in the periodic case is 

SC = M2^(6i,ci,...,6iv,c^) (16.2) 

carries two compatible local Poisson brackets, with respect to which the system BCL is 
Hamiltonian. The first ("quadratic") Poisson bracket is given by 

{bk-,bk+i]2 = —bkbk+i + Ck, {cfc,Cfc+i}2 = — CfcCfc+i 

{bk,Ck+i}2 = -bkCk+i, {cfc,&fc+2}2 = -Ckbk+2 (16.3) 

{Cfc, Ck+2}2 = —CkCk+2 
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the corresponding Hamilton function being 

Hi(6,c) 

The second (" cubic" ) Poisson bracket on BC is given by 



N 
k=l 



(16.4) 



{bk,Ck}3 = 


-Ckipkbk+i - Ck) 




{Ck,bk+l}3 = 


-Ckibkbk+i - Ck) 




{bk,bk+i}3 = 


-{bk + bk+i){bkbk+i - 


- Ck) 


{Cfc,Cfc+l}3 = 


-CkCk+iibk + bk+2) 




{bk,Ck+l}3 = 


—bkCk+i{bk + bk+i) + CfcCfc+i 


{Ck,bk+2}3 = 


—Ckbk+2{bk+i + bk+2) 


+ CkCk+\ 


{bk,bk+2}3 = 


—bkbk+ibk+2 + bkCk+i 


+ Ckbk+2 


{Ck,Ck+2}3 = 


—CkCk+2{bk+l + bk+2) 




{bk,Ck+2}3 = 


—Ck+2{bkbk+i — Ck) 




{Ck,bk+3}3 = 


—Ck{bk+2bk+3 — Ck+2) 






—Ckbk+2Ck+3 





(16.5) 



and the corresponding Hamilton function is equal to 

1 ^ 

Ho(&,c) = -^log(cfc) 

k=i 

16.2 Lax representation 

The Lax matrix of the BCL found in [ [B(J|| , |[ABC| ], is given in terms of the matrices 



(16.6) 



N N N 

U{X) = A ^ Ek+i,k, V{b, c,X) = I- bkEk,k+i + ^ CkEk,k+2 (16.7) 

k=l k=l k=l 



Theorem 16.1 The equations of motion ( |16.1| ) are equivalent to the following matrix dif- 
ferential equation: 

V = VB~AV (16.8) 

where 

N N 

A{b,c,\) = n+{UV-') = Y.bk^iEkk + \Y.Ek+i,k 

k=l k=l 

N N 

B{b,c,\) = 7c+{V-'U) = Y.bkEkk + \T.Ek+i,k 

k=l k=l 



(16.9) 
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so that also the following equation holds identically: 

UB~AU = (16.10) 

Corollary pq, [[XBq] . The matrix 

T{b,c,X) =U ■V-\b,c,X) (16.11) 
satisfies the usual Lax equation in g: 

f=[T,A] (16.12) 



This Lax equation can be given an r-matrix interpretation in the case of quadratic 
Poisson bracket, which can be lifted to an r-matrix interpretation of the Lax triads for the 
pairs {U, V) E g ® g. 

16.3 Discretization 



Since the Lax equation ( |16.12| ) has the form ( p.l|) , and moreover, the pairs (f/, V ^) satisfy 



the Lax triads of the form (|3.2| ), we can apply the recipe of Sect. ^. Taking, as usual, 
F{T) = / + /iT, we come to the discrete time matrix equation 

V = A-WB (16.13) 

with 

A = Yi+{I + hUV~^), B = Y[+{I + hV-^U) 
Moreover, since the equation 

U = A-^UB (16.14) 
holds, we have also the usual Lax equation 

f = A-^TA (16.15) 



Theorem 16.2 The equations ( p.6.13| ), ( p.6.14| ) are equivalent to the following equat 



ions: 



1 , Oik+2 , ( 1 «fc+2 \ ~ Olk+3 , . 

Ok = Ok n[ck Ck-i , Ck = Ck (16.16) 

ttfc V "fc UkUk^i ak 
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where the coefficients = ak{h, c) = 1 + 0{h) are uniquely defined by the recurrent relation 



ak = l-\ \ (16.17) 

ak-i ak-iak-2 

The following asymptotics hold: 

ak = I + hbk-i + 0{h^) (16.18) 
Remark. The auxiliary matrices A, B are bi-diagonal: 

N N 

A(6,c,A) = Y.^kEkk + h\Y,Ek+i,k (16.19) 

k=l k=l 

N N 

B{b,c,X) = Y.^k+iEkk + h\Y,Ek+i,k (16.20) 

k=l k=l 

Proof is by now standard. The general bi-diagonal structure of the factors A, B follows 
from A = + hUV^^) = n+(V^ + /if/), the latter representation implies also the recurrent 
relation for the entries of the matrix A. From (|16.14 ) one derives immediately Pk = ctfc+i- 



The equations of motion ( |16.16D follow then easily from AV = VB. ■ 

Hereafter we call the map ( p.6.16| ) dBCL. By construction, it is bi-Poisson with respect to 



the brackets ( |16.3| ) and ( |16.5| ), approximates the flow BCL due to the asymptotics ( |16.18| ), 
but is nonlocal due to the nature of the auxiliary quantities ak- 

16.4 Local equations of motion for the dBCL 

The localizing change of variables for dBCL is the map BC{h, c) i— > BC{b, c) given by the 
formulas: 

bk = hk{l + hhk-i) - hck-i, Ck = Ck{l + hhk-i) (16.21) 
As usual this is a local diffeomorphism for h small enough. 



Theorem 16.3 The change of variables (|16.21|) conjugates the map dBCL with the following 
one: 

hk{l + hhk-i) - hck^i = bfc(l + hhk+i) - hck 

(16.22) 

Cfc(l + hhk-i) = Cfc(l + hhk+2) 
Proof. Introducing the quantities 

ak = l + hhk-i (16.23) 
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we immediately see via simple check that they satisfy the recurrent relations (|16.17|) . Hence 
they represent the unique solution of these recurrencies with the asymptotics = l + 0{h). 
Now the equations of motion follow directly from (|16.16|) , ( |16.21| ), ( p.6.23| ). ■ 



Theorem 16.4 The pull back of the bracket 

{■,-}2 + h{;-h 



;i6.24) 



on BC{h, c) under the change of variables ( |16.21j ) is the following local Poisson bracket on 
i3C(b,c).- 



{bfc,Cfe} 
{bfc, bfc+i} 

{Cfc, Cfe+i} 

{bfc, c,fc+i} 

{Cfc, bfc+2} 



-/iCfc(bfcbfc+i - Cfe)(l + hhk) 
-/iCfc(bfcbfc+i - Cfe)(l + hhk+i) 
-(bfcbfe+i - Cfc)(l + hhk){l + hhk+i) 

-CfcCfc+i(l + hhk + hhk+2 + /i^(bfcbfc+i - c^) + h^ihk+ihk^, 
h{hkhk+i - Cfe))(l + hhk) 



(16.25) 



-Ck{hk+2 + /i(bfc+ibfe+2 - Cfc+i))(l + hh 



k+2 

fc+lbA;+2 ^ Cfc+1 



r/ie map ( |16.22| ) zs Poisson with respect to this bracket. 
Proof - by a straightforward but tiresome calculation. I 

Theorem 16.5 The pull-back of the flow BCL under the change of variables (|16.21|) is 

described by the following equations of motion: 



[l + hhk)ihk{h 



>fc+i — bfc_i) — Cfc + Cfc- 



;i6.26) 



Cfc = Cfc(bfc+2(1 + /ib 



fc+i) - bA:_i(l + hhk) - hCk+i + hCk-i 



Proof. We can use the Hamiltonian formalism. The pull-back we are looking for, is a 
Hamiltonian system on BC{h, c) with the Poisson bracket ( |16.25|) and the Hamilton function 
which is a pull-back of Xl^Li log(cfc). Indeed, this function is a Casimir function for 

{■, ■}2 and is the Hamilton function for BCL in the bracket h{-, -js. Calculating the equations 
of motion generated by the Hamilton function 



N 



N 



{3h)-' log(cfc) + E log(l + hhk 



k=l 



k=l 



m the bracket ( |T6l25| ), we arrive at (|16.26|) . 
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17 A perturbation of the Volterra lattice 

17.1 Equations of motion and bi— Hamiltonian structure 

Consider the following lattice system: 

ilk = Ukiwk - Wk-i + UkWk - Uk^iWk-i) 

(17.1) 

Wk = WkiUk+l - Uk + Uk+lWk+1 - UkWk) 

It may be considered as a perturbation of the Volterra lattice, therefore we adopt a name 
PVL for it. Its phase space in the periodic case is 

PV = m2^(ui,^i;i,...,mjv,w^jv) (17.2) 

The system PVL is bi-Hamiltonian. First of all, it is Hamiltonian with respect to a quadratic 
Poisson bracket on VV which is identical with the invariant quadratic Poisson bracket of the 
Volterra lattice: 

{Uk,Wk}2 = -UkWk, {Wk,Uk+l}2 = -WkUk+l (17.3) 

The corresponding Hamilton function is equal to 

N 

Y{i{u,w) = Y^{uk + Wk + UkWk) (17.4) 
k=i 

The second ("cubic") invariant Poisson bracket, compatible with the previous one, is different 
from the cubic bracket of the VL and is given by 

{Uk, Wfcjs = -UkWkiUk + Wk+ UkWk), {Wk, Mfe+lja = -WkUk+l{Wk + Uk+l) 

{Uk, Mfc+ljs = -UkUk+l{Wk + UkWk), {Wk, Wk+l}3 = -WkWk+liUk+1 + Mfe+l^fc+l) 

{Wk,Uk+2}3 = -WkUk+lWk+lUk+2 (17.5) 

The corresponding Hamilton function may be taken as 

N N 

}1q{u,w) = ^\og{uk) or Ho(m, w) = ^ log(wfc) (17.6) 

k=l fc-1 

(the difference of these two functions is a Casimir of the bracket ( |17.5| )). 
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17.2 Lax representation 

The Lax representation for PVL is given in terms of three matrices from g: 

N N 

U{u,w,X) = UkEk,k + Ek+i,k 

k=l k=l 

N 

V{u,w,\) = I - X'^^UkWkEk,k+i (17.7) 

fc=i 

N 

W{u,w,\) = I + X'^J2'^kEk,k+i 

k=l 



Theorem 17.1 The equations of motion ( 17.1 ) are equivalent to the following matrix dif- 
ferential equations: 

U = UC -AU, W = WB -CW (17.8) 
and imply also the matrix differential equation 

V = VB~AV (17.9) 

with the auxiliary matrices 

N N 

A{u,w,X) = J2i^k + Uk-iWk-i + Wk-i)Ekk + )\^Ek+i,k 

k=l k=l 

N N 

B{u,w,X) = J2(^k + UkU)k + u)k-i)Ekk + Xj2^k+i,k (17.10) 

k=l k=l 

N N 

C{u,w,X) = J2i'^k + UkWk + Wk)Ekk + Xj2Ek+i,k 

k=l k=l 

Proof - an elementary check. ■ 

It is easy to estabhsh the following fact: 

A = -K+{UWV~^), B = 'K^{V-^UW), C = -K+{WV-^U) 

so that the triples {U, V^^, W) & g ® g ® g satisfy the Lax equations of the type ( p. 21) with 
m = 3 and /(T) = T. These equations may be given an r-matrix interpretation, at least in 
the case of the quadratic bracket {■, ■}2. The corresponding quadratic bracket on g ® g ® g 
turns out to be identical with the one introduced in [|S12| . 
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17.3 Discretization 

To discretize the PVL, we can apply the recipe of Sect. |^ with F(T) = I + hT. So, we have 
to consider the following discrete time Lax representation: 

U = A-^UC, V = A-^VB, W = C'^WB (17.11) 

where 

A = n+(/ + /if/W"^), B = i[+(i + hv-^uw), c = ii+(i + hwv-^u) 



Theorem 17.2 The discrete time Lax equations ( |17.11|) are equivalent to the following equa- 
tions of motion: 

Uk = Uk—, Wk=Wk (17.12) 

where the functions ak = ak{u,w) = 1 + 0{h) are uniquely defined by the system of recurrent 
relations 

1 , Z,^ , ^ , ^(1 - h)UkWk /IV 

ak+i = I + h{uk+i + Wk) ^ (17.13) 

and the coefficients bk = bk{u,w) = 1 + 0{h), Ck = Ck{u,w) = 1 + 0{h) are given by 

ak + hukWk Uk+i-huk+i /1^71^^ 

fafc = tifc-i —r , Cfc = afc r (17.14j 

Qfc-i + hUk-iWk-i ak — hUk 

The following asymptotics hold: 

ak = 1 + h{uk + Uk-iWk-i + Wk^i) + 0{h'^) (17.15) 

bk = l + h{uk + UkWk + Wk-i) +0{h^) (17.16) 

cfc = 1 + h{uk + UkWk + Wk) + 0{h^) (17.17) 

Remark. The auxiliary matrices A, B, C are bi-diagonal: 

N N 

A = ^kEk,k + hXj2 Ek+i,k (17.18) 

k=l k=l 

N N 



B = J2bkEk,k + hXj2Ek+i,k (17.19) 

k=l k=l 
N N 

C = J2 'kEk,k + h\J2 Ek+i,k (17.20) 



k=l k=l 
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Proof. We have: A = U+{I + hUWV-^) = n+(V + hUW), since V G G_. Here 

N N N 

V + hUW = ^(1 + huk + hwk-i)Ek,k - A""^ XI ~ h)ukWkEk,k+i + hXY^ -Efc+i.fc 

k=l k=l k=l 

and the recurrent relations ( |17.13| ) for a^, the entries of the 11+ factor of this tri-diagonal 
matrix, follow immediately. The expressions for bk, ck through a^, as well as the equations of 
motion ( PTT^ ), follow directly from the Lax equations AU = UC, AV = VB, CW = WB. 



We denote the map defined in this theorem by dPVL. As usual, it shares with the system 
PVL the bi-Hamiltonian structure, the integrals of motion, and so on, but is highly nonlocal. 

17.4 Local equations of motion for dPVL 

The localizing change of variables for the map dPLV is given by the formulas 

1 + hwk-i 1 + huk nvoi^ 
Uk = Uk- , Wk = Wk- r (17.21) 

1 - /lUfc_iWfc_i 1 - ftUfcWfc 

Indeed, the following statement holds. 



Theorem 17.3 The change of variables ( |17.21| ) conjugates the map dPVL with the following 
one: 

1 + hwk-i 1 + hwk 

Ufe — = Ufc 



1 — /lUfc_iWfc_i 1 — /lUfcWfc 

(17.22) 

1 + huk 1 + huk+i 

Wfc — — ^ = Wfc 



1 - hukWk 1 - huk+iWk+i 

Proof. The statement will follow immediately, if we establish the local expressions for the 
quantities 0^: 

flfc = — 7-. 7 ^ — (17.23) 

(1 - huk-iWk-i) 



Indeed, from (|17.21|) , (|17.23|) and the formulas ( |17.14| ) we derive immediately: 

^ _ (1 + huk){l + hwk-i) 
(1 - hukWk) 



17.24) 



{1 + hUk){l + hWk) 

(1 - huk^k) 
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and then ( |17.12|) imply (|17.22|) . To prove ( |17.23D , we take this formula as a definition of 
the quantities and by means of a simple algebra verify that then the recurrent relations 
( |17.13| ) hold. The reference to the uniqueness of solution to these recurrent relations finishes 
the proof. ■ 

Theorem 17.4 The pull-back of the bracket 

{■,-}2 + M-,-}3 (17.26) 
on VV{u,w) under the change of variables ( p.7.21| ) is the following bracket on PV(u, w); 

{uk, Wfc} = -UfeWfc(l + /iUfc)(l + hwk), {wfc, Ufc+i} = -WfcUfc+i(l + hwk){l + huk+i) 

(17.27) 

The map ( |17.22| ) is Poisson with respect to the bracket ( |17.27| ). 
Proof - by a straightforward verification. ■ 

It is very interesting that the bracket ( |17.27] ) again turns out to be identical with the 
invariant local bracket ( |14.32|) of the local version of dVL (so that the contributions of 
different cubic brackets for VL and PVL are somehow compensated by different localizing 
changes of variables). 

Theorem 17.5 The pull-back of the flow PVL under the change of variables ( |17.21| ) is 
described by the following differential equations: 



life = Ufe(l + /iUfc) 



Wfc + UfcWfc _ Wfc_i + Ufc_iWfc_i 
_ 1 - /iUfcWfc 1 - /iUfc_iWfc_i 



(17.28) 



/ Ufc+i + Ufe+iWfc+l Ufc + UfcWfc \ 

Wfc = Wfc(l + /lWfe) — 

\ 1 - /lUfe+iWfe+i 1 - /lUfcWfc / 



Proof. To obtain these differential equations, one has to calculate the Hamiltonian equations 
of motion generated by the Hamilton function 

N 

(log(wfc) + log(l + huk) - log(l - hukWk) 

k=l 

with respect to the Poisson brackets (|17.27[ ). ■ 



18 Some constrained lattice KP systems 

In this section we introduce a large family of systems generalizing simulatneously the Volterra 
lattice, the relativistic Toda lattice, the Belov-Chaltikian lattice, and the perturbed Volterra 
lattice. For some reasons it is convenient to call these systems constrained lattice KP systems. 
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18.1 Equations of motion and Hamiltonian structure 

Each system of this family may be treated as consisting of m sorts of particles. The phase 
space of such systems in the periodic case is described as: 

/C^ = R"^^(i;(^),...,t;(™)) (18.1) 

where each vector 

= GM^ (18.2) 

represents the set of particles of the jth sort. We introduce the notion of the signature of 
the constrained lattice KP as an ordered m-tuple of numbers: 

e = (ei,...,ej, ej e {0,1}, ei = (18.3) 

The constrained KP lattice (hereafter cKPL(m)) with the signature e is the following system 
of differential equations: 




.0-)_„(,)IV^/„(.) _J.>, , e,(„a-4!!,)+ Y. (18.4) 

i=j + l 



Obviously, the Volterra lattice VL belongs to this class and is characterized by the sig- 
nature e = (0,0). The RTL+ flow of the relativistic Toda hierarchy also belongs to this 
class and has the signature e = (0, 1). We will see later on that BCL is Miura related to a 
simple reduction of cKPL(3) with e = (0, 1, 1), while PVL itself is a less trivial reduction 
of cKPL(3) with e = (0, 1,0). The system cKPL(m) with the signature consisting only of 
zeros, e = (0, . . . , 0) is nothing but the Bogoyavlensky lattice BLl(m — 1). Also the v^}^ = 
reduction of cKPL(m) with the signature e = (0, 1, 1, . . . , 1) coincides with BLl(m — 1). 

Theorem 18.1 The system ( |18.4|) is Hamiltonian with respect to the following Poisson 
bracket on /Cm-' 

{4^^4^i}2 = -9t;?4^i (18.5) 
b^.^\ = -^fv\:\ {4^\4'ii}2 = -4^Mti for l<^<J<m (18.6) 
with the Hamilton function 

m N 

Mv) = Hi(t;«, . . . , .;(-)) = J2i: 4^'^ (18.7) 

j=i k=i 

(independent of the signature). 

This statement can be easily checked and generalizes the quadratic brackets ( |14.4| ) and 
( |15.5| ) for the VL and the RTL, respectively. It would be important to find out, when do 
the analogs of the linear bracket (for the RTL) and of the cubic bracket (for both the VL 
and the RTL) hold, and to find the corresponding expressions. 
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18.2 Lax representation 



The natural Lax representation of the cKPL(m) ( |18.4D hves in and is given in terms of 
the following matrices from g: 

N N 

f/i(t;«,A) = \Y.Ek+i,k + T.^k'^Ek,k (18.8) 

fc=l k=l 
N 

Vj{v^^\X) = I + cjj\-'Y.^k^Ek,k+i, j = l,2,...,m (18.9) 

k=l 

where 

a, = |_| ^J:M = l-2e, (18.10) 



Theorem 18.2 The equations of motion ( |18.4| ) are equivalent to the following matrix dif- 
ferential equations: 

ill = UiBm-BiUi (18.11) 

r V,B,_, - B,V, e, = 0] 

Vj = I \ 2 < j <m (18.12) 

where 

N j m N 

B,{v,X) = T.{T.^k^+ E vti)Ek,, + XT.Ek+i,k (18.13) 

k=l i=l i=j+l k=l 

The evolution of the monodromy matrices 

Tj{v^'\ . . . , A) = Vp{X) ■ . . . ■ VpiX) ■ U,{X) ■ ^^-(A) ■ . . . ■ V;^rW (18.14) 
is governed by usual Lax equations: 

Tj = [T„B,] (18.15) 
The matrices Bj allow the representation 

= n+{Tj) (18.16) 
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Proof - an easy check. ■ 

It is easy to see that this Lax representation is exactly of the form (|3.2| ), if one considers 
the m-tuple of matrices 

as the Lax matrix. This Lax representation allows an r-matrix interpretation for the 
quadratic bracket of the previous theorem. As a matter of fact, the corresponding quadratic 
bracket in gr®™ literally coincides with the bracket introduced for the BLl in [ S12|| . The Lax 
representation (|18.11|) , ( |18.12|) serves as a starting point for applying the recipe of Sect. |. 



18.3 Discretization 

Taking in this recipe F{T) = / + hT, we come to the discrete time Lax equations 




e, =0 



6, = 1 



2 < j < m 



[18.17) 



with Bj =U+{I + hTj). 

Theorem 18.3 The discrete time Lax equations ( |18.171 ) are equivalent to the map v \—>- v 
described by the equations 



Vk -^k (1) 



V 



.0') 



V 



6, =0 



e, = l 



(2 < J < m) 



(18.18) 



where the functions b^^^ = b^k\v) = 1 + 0{h) satisfy the following equations: 



.0-1) 



''fc+1 ''-^fc+l 



(1) 



,(1) 



(1) 



bk' - 
Aj) _ 1^ U) 



4" 



bl"" — hv 



k 



,0-1) 



0) 

k 



0, 2 < j < m 



e,- = 1, 2 < j < m 



[18. 19) 
[18.20) 
[18.21) 
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Proof. First of all notice that the matrices Bj must have the following structure: 

N N 
k=l k=l 

Now the equations of motion are derived straightforwardly. For example, the = 1 variant 
of the last m — 1 equations in ( |18.17| ), i.e. the matrix equation Bj^iVj = VjBj, is equivalent 
to the following system of scalar equations: 



This is equivalent to the corresponding variant of equations of motion (|18.18|) together with 
the relation ( |T8:21| ). ■ 



Remark. It is important to notice that the statement of the last theorem deviates 
from the usual scheme in that it does not contain a system of equations which determine 
b^'''* uniquely. In fact, in order to find such a system one has to deduce from the equations 
( |m9| )- ([lMT|) m formulas of the type 



^k+q _ ^k 



Here the number q does not depend on j; ^'^■'^ are certain expressions of the type Jli^i '4'k+nii 

and all iIj^^^ = b^*'' — hv^^^ or ijj^^^ = (^b^* + hv^^^^ . This, in turn, implies that there hold 
certain equations of the type 

b^^ • . . . ■ bg,_i*l-'') = const (18.22) 

(here we assumed for definiteness that g > 0). The value of the constant on the right-hand 
side is uniquely defined by the conditions Bj = + hTj). (As a matter of fact, it is easy 
to see that this constant does not depend on j). The value of the constant being determined, 
the equations ( p.8.22| ) give the desired system which defines b^j^^ uniquely. However, the outfit 
of this system depends heavily on the signature e of the cKPL, and the general formulas 
would contain too many indices to be instructive enough. It is simpler to derive such formulas 
for each concrete signature separately. Nevertheless, the formulas (|18.19|) -( |18.21j ) completely 
characterize the coefficients of the matrices which serve as the factors Ilj^{aI+hTj) with some 
a, so that every solution of this system leads to a discretization based on the factorization of 
I + h'T with h' = h/a, which enjoys all the positive properties of our general construction. 
We call the maps introduced in the previous theorem dcKPL(m). 
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18.4 Local equations of motion for dcKPL 

The dcKPL(m) can be brought into the local form for an arbitrary signature e. 
Theorem 18.4 The change of variables ICmiy) ^ KLmiy), 

i-i 



V 



V? n (l + ^v«) ■ (l + -11(1 + i) (18-23) 



i=l 



i=j+l 



conjugates dcKPL(m) with the following map: 

i-i 



V? n (i + ^v«) • (i + e.h^E) ■ n (i + ^vEi 

i=l *=i+i 

=vFn(i+^v«,).(i+6,/.vg,). n (i+/^v« 



i=j+i 



(18.24) 



Proof. It is easy to calculate that if ( |18.23| ) holds, and if the quantities b^"''' are defined by 
the formula 

j 



bi^-) = n (i + /^v«) n (i+/^v«i 



^18.25) 



1=1 



i=j+i 



then 



i=l 



+ z^^^? = n (i + h^f) n (i + ^vEi 



j=i 



6,=0 



and it is easy to check now that the equations ( |18.19D -( p.8.21| ) are satisfied. Indeed, for = 
we find: 



- hv\:^ bg, - h^i 



1 + hv^i^ 



which proves (|18.20|) , while for = 1 we find 



(i) 



which proves (|18.21|) . The verification of ( |18.19| ) is completely analogous. The pull-back of 
the equations of motion is calculated now straightforwardly. ■ 
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Unfortunately, we do not now a general formula for the second invariant Poisson structure 
for cKPL's. This prevents us from applying our general scheme for finding the local invariant 
Poisson brackets for the localized maps. However, by a direct analysis of equations of motion 
the following statement can be proved. 
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Theorem 18.5 The pull-hack of the system ( p.8.4D under the change of variables ( |18.23D 
given by the formula 

Cj-l m 
n (i + ^viti) ■ (i + e./^vgi) ■ n (i + 
i=l «=i+i 

i— 1 m \ 

- n (l + h^f) ■ (l + ^Mf-i) -11(1 + ^vEi) /h (18.26) 

1=1 i=i+l / 

18.5 Example 1: e= (0,0,0) 

As illustrations we consider the systems with m = 3 - the simplest possible ones after VL 
and RTL+. We denote for simplicity 

(1) (2) (3) 

The Hamilton function is in all cases equal to 

N 

Hi(M,t;,w) = ^{uk + Vk + Wk) 

k=l 

The nonvanishing Poisson brackets consist of the signature independent part, 

{Uk,Vk}2 = -UkVk, {Vk,Uk+l}2 = -Uk+lVk 

{Uk,Wk}2 = -UkWk, {Wk,Uk+l}2 = -Uk+lWk (18.27) 

{Vk,Wk}2 = -VkWk, {Wk,Vk+l}2 = -Vk+lWk 

supplemented by 

for those j where = 1. 

In particular, if all = 0, then all nonvanishing Poisson brackets of the coordinate 
functions are exhausted by ( |18.27D , and we arrive at the system 



Uk = Uk{Vk + Wk - Vk-l - Wk-l) 

Vk = Vk{uk+i + Wk-Uk- Wk-i) (18.28) 

Wk = WkiUk+l + ffc+l -Uk-Vk) 
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which becomes the usual Bogoyavlensky lattice BL1(2) after the re-naming 



The localizing change of variables for its discretization: 



Uk = Ufe(l + /lVfc_i)(l + /lWfc_i) 

Vk = Vfc(l + /iUfc)(l + /iWfc_i) (18.29) 
Wk = Wk{l + huk){l + hvk) 

The local discretization of the system ( |18.28D : 

Ufe(l + hYk-i){^ + hwk-i) = Ufc(l + /iVfc)(l + hwk) 

Vfc(l + huk){l + hwk-i) = Vfc(l + /iufc+i)(l + hwk) (18.30) 

Wa.(1 + /iUfc)(l + hvk) = Wfc(l + /iUfe+i)(l + hvk+i) 

18.6 Example 2: e = (0, 1, 0) 

In this example the signature dependent part of the Poisson brackets reads 

{Vk,Vk+l}2 = -VkVk+1 

and the equations of motion take the following form: 

Uk = UkiVk + Wk- Vk-l - Wk-i) 

Vk = Vk{uk+i + Vk+i + Wk - Uk - Vk^i - Wk-i) (18.31) 

Wk = Wk{Uk+l + Vk+l -Uk -Vk) 

The localizing change of variables for the discretization of this system: 

Uk = Ufc(l + /lVfc_i)(l + /lWfc_i) 

Vk = Vfc(l + /iUfc)(l + /iVfc_i)(l + /iWfc_i) (18.32) 
Wk = Wk{l + huk){l + hvk) (18.33) 



The local form of equations of motion for the discretization of ( 18.31|) : 



Ufc(l + /iVfc_i)(l + hwk-i) = Ufc(l + /iVfc)(l + hwk) 

Vfc(l + /iUfc)(l + /iVfc_i)(l + hwk^i) = Vfc(l + /iUfc+i)(l + /iVfc+i)(l + hwk) 
Wfc(l + huk)il + hvk) = Wk{l + huk+i){l + hvk+i) (18.34) 
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Let us mention that the system ( |18.31| ) allows an interesting reduction 

Vk=UkWk (18.35) 

which is, moreover, compatible with the quadratic Poisson brackets. In this reduction we 
arrive at the system PVL. It is easy to check that in the variables u^, Vfc,Wfc the reduction 
( |18.35| ) takes the form = UfcWfc(l + hv^), so that 

UfcWfc 1 

Vfe = , 1 + 



1 - /lUfcWfc 1 - /lUfeWfe 



This makes a link with the results of Sect. 17 



18.7 Example 3: e = (0, 1, 1) 

In this case the Poisson brackets ( |18.27D have to be supplemented by 

{ffe, ffc+i}2 = -VkVk+1, {Wk, Wk+l}2 = -WkWk+l 

and the equations of motion take the form: 

ilk = Uk{vk + Wk - Vk-i - Wk-i) 

Vk = Vk{uk+i + Vk+i + Wk - Uk - Vk-i - Wk-i) (18.36) 

Wk = Wk{Uk+l + Vk+l + Wk+l - Uk - Vk - Wk-l) 

The localizing change of variables for the discretization of this system: 

Uk = Ufc(l + /lVfc_i)(l + /lWfc_i) 

Vk = ^rk{l + huk){l + hvk-l){l + hwk-l) (18.37) 
Wk = Wk{l + huk){l + hvk){^ + hwk-i) 



The local form of equations of motion for the discretization of (|18.36|) : 



Ufc(l + /iVfe_i)(l + hwk-i) = Ufc(l + hvk){l + hwk) 

Vfc(l + /iUfe)(l + /iVfc_i)(l + hwk-i) = Vfc(l + /iUfc+i)(l + /iVfc+i)(l + hwk) 
Wfc(l + huk){l + hvk)il + hwk-i) = Wfc(l + /iUfc+i)(l + /iVfe+i)(l + hwk+i) 

(18.38) 



Let us discuss the following reduction of the system (|18.36|) : 



Uk = (18.39) 
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It is compatible with the quadratic Poisson brackets, so that we arrive at the following 
reduced system: 

Vk = Vk{Vk+l + Wk- Vk-l - Wk-l) 

(18.40) 

Wk = WkiVk+l + Wk+l - Vk - Wk-l) 

Interestingly enough, this system is again nothing but the usual Bogoyavlensky lattice 
BL1(2), which becomes obvious after the re-naming 

Vk ^ a2k-i, Wk ^ a2k 

So, we have found the third Lax representation for BL1(2). 

It is easy to see that the maps A4± : JC^^O, v, w) BC{h, c) defined as 

M+ : h = Vk + Wk, Ck = VkWk+i (18.41) 

and 

M- : bk = Vk+i + Wk, Ck = Vk+2Wk (18.42) 



conjugate the flow ( 18.40 ) with the Belov-Chaltikian flow BCL, and are Poisson, if the both 
spaces are equipped with the brackets {■, ■}2- So, the system BCL is Miura related to BL1(2) 
(this fact is similar to the Miura relation between the Toda and the Volterra hierarchy). 

The localizing change of variables for the discretization of the reduced system ( |18.40| ) is 
given by 

Vk = Vk{l + h\k-i){l + hwk-i), Wk = Wk{l + h\k){l + hwk-i) (18.43) 
and the corresponding local equations of motion read: 

Vfc(l + hvk-i)il + hwk-i) = Va.(1 + /iVfc+i)(l + hwk) 

Wfc(l + /iVfc)(l + hwk-i) = Wfc(l + /iVfc+i)(l + hwk+i) 



;i8.44) 



So, the discretizations of BL1(2) based on different Lax representations, agree with one 
another. 

It turns out that the Miura maps are still given by nice local formulas, when trans- 
lated to the localizing variables. Namely, the following diagram is commutative: 
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M 



± 



/C3(0,V,W) 



BC{h, c) 



}Cs{0,v,w) 



M- 



i3C(6, c) 



if the maps M± are defined by the formulas 

M+: l + /ibfc = (l + /ivfc)(l + /iwfc), Ck = VkWk+i{l + hwk){l + hvk+i) (18.45) 
and 

M_ : 1 + /ibfc = (1 + /iVfc+i)(l + /iWfc), Cfc = Vfc+2Wfc(l + /iVfc+i)(l + /iWfc+i) (18.46) 
This statement may be verified by a simple calculation. 



19 Bruschi— Ragnisco lattice 

The Bruschi-Ragnisco lattice (hereafter BRL) was introduced in ||BR2 

h = h+iCk - hck-i, Ck = Ck{ck - Cfc„i) 



:i9.ii 



It may be considered either under open-end boundary conditions (^at+i = Cq = cat = 0), or 
under periodic ones (all the subscripts are taken (mod A^), so that Cq = Cjv, ^at+i = bi). The 
phase space of the Bruschi-Ragnisco lattice: 



(6l,Ci, . . . ,bN,CN) 



Two compatible brackets may be defined on BTZ such that the system ( |19.1|) is Hamiltonian 
with respect to each one of them. The linear Poisson bracket is given by 



{bk,Ck}o — — {^fe+i,Cfe}o — —Ck 



:i9.2) 
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while the quadratic one - by 

{bk,bk+i}i = -h+iCk, {bk,Ck}i = cl, {bk,Ck+i}i = -CkCk+i (19.3) 
The corresponding Hamilton functions are: 

N N 

Hi (6, c) = ^ bk+iCk and Ho(6, c) = ^ 6^ (19.4) 

fc=l k=l 

for the brackets {■, - jo and {■, -ji, respectively. 

It has been pointed out in |PR2|| that this system allows a very complete study by 
different methods of the soliton theory. As was demonstrated in [ p^ , this is due to its 
extreme simplicity. Namely, in a certain gauge the Lax representation of this system is a 
linear matrix equation. Namely, if the entries of the Lax matrix T = T(6, c) G gr = gl{N) 
are defined as 

bj UlZl Ci k<j 
Tkj = { ^ (19.5) 

bj (nt-^ k > J 



then the system ( |19.1|) is equivalent to the equation 

T = [T, M] (19.6) 

with the constant matrix M 

N-l N-1 

M=Y^ Ek,k+i or Ek,k+i + CEn^ (19.7) 

k=l k=l 

for the open-end and periodic case, respectively (in the latter case it is supposed that the 
dynamics of the BRL is restricted to the set Ci . . . = C) . The whole hierarchy of the BRL 
consists of equations 

T = [T, M™] 

which are linear and may be immediately integrated: 

Tit) = exp{-tAT") ■ T(0) • exp(tM'") (19.8) 

The brackets ( 19.2|) , (|19.3|) were shown in |^2|| to give a coordinate representation of certain 
Lie-Poisson brackets on g, restricted to the subset of the Lax matrices T{b, c). 

Obviously, the recipe of Sect. H cannot be literally applied to the BRL. However, the 
philosophy behind this recipe is, of course, apphcable, and requires to seek for the discrete 
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time Bruschi-Ragnisco lattice in the same hierarchy. It should share the Lax matrix with 
the continuous time system, and its explicit solution should be given by 

T(n/i) = (/ + /iM)-"T(0)(/ + /iM)" (19.9) 

(cf. (|19.8| )). Hence the corresponding discrete Lax equation should have the form 

f = {I + hM)-^T{I + hM) (19.10) 



Theorem 19.1 The discrete time Lax equation ( |19.10 ) is equivalent to the following map 
on the space BTZ: 

hk{l + hck-i) =bk + hbk+iCk, 4 = Ck ^^!^~^ (19.11) 

i + riCk-i 

Proof - an easy calculation. ■ 

By construction, this map is Poisson with respect to the both brackets ( |19.2|) , (|19.3|) . 
We see that the extreme simplicity of the BRL allows to find its local discretization in the 
original variables. The localizing change of variables is not necessary for this system. 



20 Conclusion 

This paper contains a rich collection of examples illustrating the procedure of constructing 
local integrable discretizations for integrable lattice systems. The construction is based on 
the notion of the r-matrix hierarchy and consists of three steps of a rather different nature. 

The first step is to find a Lax representation for a given lattice system, living in an 
associative algebra g. This Lax representation has to be a member of a hierarchy governed 
by an /^-operator on g satisfying the modified Yang-Baxter equation. In all examples treated 
here this operator is simply a difference of projections to two complementary subalgebras. 

The second step is an application of a general recipe for integrable discretization. This 
step is almost algorithmic, the only non-formalized (and, probably, non-formalizable) point 
being the choice of the function F{T) approximating exp(/iT) for T & g (cf. Sect. 
In all examples treated here the simplest possible choice F{T) = I + hT works perfectly. 
The difference equations obtained on this step share the invariant Poisson structures, the 
integrals of motion, the Lax matrices, etc. with the underlying continuous time systems. 
However, as a rule, they are non-local. This feature is unpleasant from both the esthetical 
and the practical point of view, because it makes the equations ugly and not well suited for 
practical realization. 
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The third step is finding the locahzing change of variables. This step is again absolutely 
non-algorithmic (at least, at our present level of knowledge). These changes of variables 
have remarkable properties: they often produce one-parameter local deformations of Poisson 
brackets algebras, and always produce one-parameter integrable deformations of the lattice 
systems themselves. At the moment we cannot provide a rational explanation neither for 
these properties nor for the mere existence of localizing changes of variables. However, our 
collection seems to be representative enough to convince that these phenomena are very 
general. We feel that they are connected with the Poisson geometry of certain r-matrix 
brackets on associative algebras and of monodromy maps, but we prefer to stop at this 
point. 
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